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ABSTRACT 

The parameters of the) efflux from a helium dewar in space were 
numerically calculated. The flow was modeled as a one-d1mensional , 
compressible, ideal gas with variable properties. The primary boundary 
conditions are flow with friction and flow with heat transfer and friction. 
Two PASCAL programs were developed to calculate the efflux parameters: EFFLUXD 
and EFFLUXM. EFFLUXD calculates the minimum mass flow for the given shield 
temperatures and shield heat inputs. It then calculates the pipe lengths, 
diameter, and fluid parameters which satisfy all boundary conditions. Since 
the diameter returned by EFFLUXD is only rarely of nominal size, EFFLUXM 
calculates the mass flow and shield heat exchange for given pipe lengths, 
diameter, and shield temperatures. 
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NOMENCLATURE 


A 

B 

S 

D 

f 

h 

•‘f 

L 

M 

m 

Nu 

P 

Po 

Pp 

Qw 

q 

Rej- 

St 

T 

I'd 

Tw 


aw 


V 

X 


area 

function of M (* [1 + (y -l)/2] M^) 

specific heat 

diameter 

friction factor 

fluid film constant 

fluid thermal conductivity 

pipe length 

Mach number 

mass flow 

Nusselt number 

pressure 

stagnation pressure 
Prandtl number 
shield heat exchange 
heat flux 
Reynolds number 
Stanton number 
temperature 
stagnation temperature 
wall (shield) temperature 
adiabatic wall temperature 
velocity 
distance 


Greek Letters 


Y 

e 

P 

P, 


specific heat ratio 
shield effectiveness 
density 

stagnation density 
viscosity 


Subscripts 

1 entrance 

2 exit 
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1. INTRODUCTION 

One of the primary objectives in the design of a long lifetime helium 
dewar is to minimize the parasitic heat load. This is accomplished to a large 
extent, through low-conduction supports and vapor-cooled shields imbedded in 
multi-layered jnsulati on (MLI). The vapor-cooled shields can be modeled as a 
series of heat exchangers with constant temperature walls and constant heat 
input to the helium efflux. With the ambient pressure of space being 
approximately zero, the fluid flow is compressible with a Mach number of unity 
at the exit throat. The solution of this model basically consists of 
determining (1) the length of pipe required to provide the necessary heat 
transfer and (2) the diameter required for the optimum mass flow to reach the 
Mach number of unity at the throat. 

This solution entailed the development of two PASCAL programs: EFFLUXM 
and EFFLUXD. EFFLUXD is used to optimize mass flow and define the ler?gth and 
diameter of the heat exchanger pipe. Since this diameter is only rarely and 
by chance of nominal size, EFFLUXM is used to determine the mass flow and heat 
absorption of a similar heat exchanger but with a nominal diameter close to 
theoretical diameter determined through EFFLUXD. These programs ara explained 
in APPENDIX A. 

Although the original scope of the project was to develop these programs 
with the number of shields, the temperature of these shields and their 
position as given, a second-law thermal analysis of the shield system would 
seem beneficial to the project as a whole. Preliminary calculations indicate 
that the optimal number of shields is '(ess than the five previously cited by 
other sources. Information on this subject will be forwarded at a later 
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2. THEORETICAL CONSIDERATIONS 

The shield heat exchanger system is a bit unusual in that both shield 
temperature and heat input are specified at each shield from an overall heat 
transfer analysis of the entire insulation system developed elsewhere. If a 
shield heat exchanger system is analyzed with the effectiveness method, it 
becomes clear that for a given effectiveness the mass flow is a dependent 
variable. This mass flow may be minimized as an explicit function v.f shield 
temperature, heat input, and effectivaness. As an example of the method, the 
three-shield system shown in Fig. 1 will be analyzed. The initial stagnation 
temperature, the shield temperatures, the shield heat inputs, and the shield 
effectiveness (or shield-to-hel ium exit temperature difference) are all 
given. From the definition of effectiveness as shown in Eq. (1), the fluid 
stagnation temperatures at points 1, 2, and 3 may be calculated. 

® “ ^^0 out ” ^0 in^'^^^w “ ^0 in^’ 

Note that the stagnation temperature at points 1 and 2 are outlet conditions 
for one shield and inlet conditions for the next. A single mass flow, 
however, will generally create different effectiveness in each shield. 
Consequently, the optimal mass flow must be found which will provide the 
maximum effectiveness in one shield. To accomplish this objective, each jth 
mass flow is calculated as a function of the initial stagnation temperature, 
the jth wall temperature, and the summation of heat inputs to the jth wall. 
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A, = [r q^J/[c„(T„, - T 
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The maximum resulting mass flow is the optimal mass flow .for the given 
conditions. With this mass flow, the actual shield effectivenesses can be 

calculated. For constant fluid properties, these effectivenesses would 
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Figure 1 Temperature profile of helium pipe flow 
over three radiation shields. 
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directly define the heat exchanger lengths. The widely varying properties of 
the helium flow, however, force the sectional solution used in the EFFLUXD 
program. 

In addition to the problem of variable fluid properties, the flow Is 
compressible with a Mach number of 1 at the exit throat. Since the mass flow 
and length have already been set, the cohyition of M » 1 at the exit must be 
satisfied by determining the proper diameter. This Is done implicitly in 
EFFLUXD as shown in APPENDIX A. If a nominal rather than calculated diameter 
Is used, the mass flow Is no longer optirnized and must be varied until the 
exit condition of M * 1 is reached. This is done implicitly in EFFLUXM as 
shown in APPENDIX A. 

Both of these programs use a number of functions and procedures which 
solve for the exit Mach number of a short section of pipe, given the inlet 
Mach number and the boundary conditions. The first of these procedures 
calculates the pipe entrance Mach number as a function of mass flow and 
diameter as explained in APPENDIX B. 

For sections with friction only (connecting pipes between shields or 
valves), the flow is assumed laminar below a Reynolds number of 2000 and 
turbulent above that. This transition changes the frictional and thermal 

behavior of the fluid. In laminar flow, the friction factor is a function of 
Reynolds number only and the Nusselt number is constant; in turbulent flow, 
both the friction factor and Nusselt number are functions of Reynolds and 
Prandtl numbers as illustrated more clearly in APPENDIX C. 

Frictional effects are also important in a heat transfer section of 
sufficient length. Before the exit Mach number of a short section may be 
calculated, the exit stagnation temperature must be known. This can be 
calculated with the simple energy balance in APPENDIX D. One interesting 
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aspect of this balance is that it is a function of the adiabatic wall 
temperature. For laminar flow, this is equivalent to the stagnation 
temperature; for turbulent flow, the adiabatic wall temperature lies between 
the fluid stagnation temperature and the actual fluid temperature. However, 
for subsonic pipe flow, the adiabatic wall temperature can be assumed to equal 
the stagnatitm temperature with no appreciable error. 

After the exit stagnation temperature is calculated, the governing 
differential equation for compressible flow with heat transfer and friction 
may be simplified using Reynolds analogy and the previous energy balance as 
described in APPENDIX E. The resulting equation is solved as shown in 

APPENDIX F with a fourth-order Runge-Kutta method. 

Within each small section, the fluid properties may be assumed constant, 
but overall 'thermal conductivity, viscosity, and Prandtl number vary as a 
function of temperature. The Prandtl number, as a function of absolute 

temperature, was modeled with a fifth-order least squares fit as shown in 
Fig. 2, while the thermal conductivity and viscosity are modeled with third- 
order fits as shown in Figs. 3 and 4. The temperatures used are the section 
inlet fluid temperatures. 

The diameter, once determined, remains constant. Any constriction of the 
diameter in beginning sections increases the diameter in later sections to 
maintain the optimal mass flow. Furthermore, diameter increases would either 
require the design and manufacture of smooth venturi sections or would result 
in pressure loss across the area change. This added pressure loss would again 
require increasing the diameter even more or require increasing the mass 
flow. Earlier suggestions of diameter increases seemed to arise from the 
problem of extremely long pipes; heat transfer sections about two orders of 
magnitude longer than needed were proposed since, at the time, the heat 


FIGURE 2 



Figure 2 Prandtl Number versus Temperature 





Figure 3 Thermal Conductivity 
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Figure 4 Viscosity versus Temperature 
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transfer mechanism was not fully understood. Even though It is difficult to 
see the advantage of area changes, a procedure to deal with them is in EFFLUXD 
as shown in APPENDIX H. This approach, however, was never utilized. 

To’' summarize the design procedures, the following components can be 
identified: 

1. The sections in contact with heat shields are treated as heat 
exchangers with both heat transfer and friction influencing the flow 
of gas. 

2. The connecting sections between heat shields are considered as 
adiabatic lengths with friction only. 

3. Valves in the system can be treated as adiabatic sections as in point 
2, with equivalent lengths adjusted to yield the proper pressure drop 
for a given flow rate and entrance conditions. Depending on the data 
available for the valve, this calculation will require iterations by 
the user of the program. 
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APPENDIX A 

The Inputs and outputs of the PASCAL programs ^FFLUXD and EFFLUIW should 
be self-explanatory. The units are consistently SI. Example runs of both 
programs’’ with s five-shield system follows. The input data used were taken 
from "Feasibility Study for Long Liftime Helium Dewar," Parmley, R. T. , 
Lockheed Missiles and Space Company, Inc., Contract No. NAS2-10848, NASA, 


ORIGINAL PAGE IS 
OF POOR QUALITY, 


Hfi.cmuio 

HOCiiH crrLuiD( mm ,ovTrirr,irroaTi).trrnot .crrimi .rovro) ; 





•it 




ttt 

(tn 

tn 

tit 

tti 

(tti 

tn 

ttt 

tn 

<tit 

tn 

ttt 

tn 


in 

in 

III 

(ttt 

ttt 

III 

III 

(til 

ttt 

III 

III 





(iiiittitt 



tn 




III 

(tit 

in 

III 

in 

(tn 

tn 

III 

III 

(ttt 

tn 

in 

III 

(ttt 


tn 

III 

(ttt 

ttt 

ttt 

tn 


ttt 

ttt 

tti 


ttt 

ttt 

nil 


ttt 

tn 

tttl 


III 

tn 

ttt 

III 


I 

tn 

in 

ni 

tn 



III 

III 

III 

III 


tt ) 

tn 

ni 

III 

III 

III 

till 

III 

III 

III 

III 

III 

till 

III 

III 

III 

III 

tn 

III} 

III 

III 

III 

III 

III 

in| 

tn 

tn 

III 

III 

III 

III) 

III 

tti 

III 

III 

in 

III) 

tn 

III 

III 

III 

in| 

tn 

III 

III 

III 

III) 

ttt 

III 

III 

III 

III 

III) 

III 

III 

III 

III 

III 

III) 

III 

m 

III 

III 

III 

III) 

tn 

tn 

III 

III 

III 

in| 

III 

tn 

in 

ttt 

III 

III) 

III 

in 

tti 

III 

III 

III) 

tittiti 

ttiiii 

ttt 

tn 



nini 

itttt 

tn 

in 



ttni 

tttl 

ttt 

III 




(titttttttttf tttttttttltftlttttlttttltttflllltttltltttttttttfttttlltttMttttlll 

(tfltttttttfttfltxtftifttttttttttKftttti j. c. CKATO IND I. W. IRINDUT 
(tffiitittttttttttttttcttatttttfttf nut 
(ttittttttttnttttctiiitntttttttttnttt 

(ntnitinnntttntntttntttnntttn WII?EISin OF ItllHOIS 

(ntttitnttttttttiintnttttttntntttt umeI COMTXICT OF mi, IMES lESEliCE 
(tttnntttttttttifitttttttttnttttittti 

(ttntnfitttnntttftntnttttttnttttt JUIT iti: 

(•tnttnttntntftintnitfnttttitttitntttattittntttintttittttntttnttni 

(ttnntittttnttnnttnititttitntntttttitiKtntninititntttttntntttttti 




(* THIS PROCRAH DETERNIHES CERTAIH FARAHETERS OF THE EFFLUX RESULTING*) 
(• FIOR THE PARASITIC HEAT LOAD ON A HELIUM OEVAR IN SPACE. THE SHIELD *) 
(• STSTEK CONSISTS OF CONSTANT 'fEHPERATUHE SHIELDS VITH CONSTANT HEAT *) 
(• INPUTS AS PREDETERMINED BT THE THERMAL CHARACTERISTICS OF THE *) 
(• SATELLITE. THE MINIMUM MASS FLOW WHICH CAN SATISFY THESE INPUTS IS *) 
(* FIRST CALCULATED. THIS MASS FLOW IMPLIES THE LENGTH OF THE SHIELD *) 
(* HEAT EXCHANGERS. THE MASS FLOV, PIPELENTHS, AND THE NECESSARY *) 
(• CONDITION OF CHOKED FLOV AT THE EXIT ARE THEN USED TO CALCULATE THE *) 
<■ PIPE OIANETEH. THIS PROCRAM HAY BE USED IN CONJUMCTIOKvVITH EFFLUXM *) 
(• WHICH SOLVES FOR MASS FLOV GIVEN PIPE LENGTHS AND A NOMINAL DIAMETER *) 


CONST 

PM.H1SPIAS3; 

l>»77:(*JfKC-I«) 

I>1.47: 

CPtSXBI; 

nP8 

DIVIS10NS*ARRAYCI.. 20001 OF REAL; 
COUKT>ARRATC0..201 OF REAL; 
ICODHTsARRAYC&..20I OF INTEGER; 
REALOUTaFILE OF REAL; 

INTOUTaFIlE OF INTEGER; 


NUMVALLS , NUMDI V , SFACECKT ,X ,MTEST: INTEGER ; 

I .lOOPCOUNF , PREY , TAB . P IPECXT, VALLCNTNTECER; 
EFFOUTD.POUTD;TEn; 

EFFIN01:£EAL0UT: 

EFFINDX:INT0UT; 

PJ4HO.DIA:OIVISIONS; 

X,OX,H,PO,TO.Q,RE,F,NV;OI7ISIOKS: 

OV, TV, LTEMP . LSPACE , EFF , OVALL : COUNT ; 
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1. ITOm . lOVIOUKO , MI BOUND . T8SAVC ; RtU ; 

n>OT,IUXErr.OUI,OIVlENCTM:IEll; 

Ht.Tat,rBl,Tl.ri,IHOi.ri.lEt,Ql,aXl,NUI;REAli 


PIOCEOOIE 1NDATA(INDEZ:IKTECEX): 

BEGIN <• IKTEBACTIVE INBUT <) 
ir INDEX ) I THIN 
BEGIN 

(t BEAD IN OITA •) 

niTElNC INPUT DATA AS PIOHPTEO BT nBMINAL'); 

VliniNC '); 

niTElNCVHEN ASEEO POM A SERIES OF DATA SUCH AS THE SHIELOM; 
UBITEINCTEKPEIATURES. INPUT THE DATA IN THE ORDER FROM DEVU'); 
VIITEINCTO THE OUTER BOUKDART.Mj 
VRITEINC M; 

VRITEIKC ENTER THE NUKIER OF RADIATION SHIELDS*): 

REAOIH; 

RCAO(NUXVAllS); 

VHITE(EFFIND2,NUKVALIS); 

VRITELNI 'NUMBER OF SHIElOSi'.HURVAllS:S); 

VRITELNC '); 

VRITEINI 'ENTER SINGLE DIVISION LENGTH (N)'); 

REAOLN: 

READ(OITLENGTM) ; 

VRITEICFFIKDl.OITLENCTH); 

VRITELNC 'DIVISION LENGTH (N)-'.0IVLXNCTH:1S:S); 

VIITELHC 

VRITELNC 'ENnX PIPE LENGTHS NOT ON SHIELDS CK)'); 

FOl I:al TO MUHVALLSfl DO 
lECIH 
lEAOLH: 

REAOCLSPACCm); 

VRITECEFFINDI.LSPACZII]); 

END:C<rOR*) 

VRITELNC 'ENTER TEMPERATURES OF UOIATION SHIELDS CE)‘); 

FOl I:>1 TO HUKVALLS DO 
lEClM 
REAOLN: 

READCTUCIl); 

VRITECEFFINDI.TVni): 

END;C< FOR *) 

VRITELNC 'ENTER HEAT INPUTS OF RADIATION SHIELDS CVI'I; 

FOR I:>l TO KVKVALLS 00 
BEGIN 
REAOLN; 

READCOVIIl); 

VHITECEFFINDl.aVdl); 

END: C 'FOB*) 

VRITELNC 'ENTER HAZIHUM HEAT EXCHANGER EFFECTIVENESS AS A FRACTION'): 
BEADLN; 

REAOCHAXEFF); 

VHITECCFFtNOMUXEFF): 

VRITELNC 'SHIELD EFFECTIVENESS* ' .MAXEFF; IS: 7); 

VRITELNC '); 

VRITELNC 'ENTER INITIAL STAGNATION T CX)'); 

EEAOLN; 

AEADCTOm): 

VRITECEFFINDUBII]): 

TV(Bl:xTOm; 

VRITELNC 'INITIAL STAGNATION TEMPERATURE CI)t ' jgni:lS:S): 

VRITELNC '); 

VRITELNC 'ENTER INITIAL STAGNATION P CPA)'); 

REAOLN: 

AEADCFOni): 

VRITCCEFFINDhPlII]); 

VRITELNC 'INITIAL STAGNATION PRESSURE CPA)*'.P0C11:1S:S); 

VRITELNC '); 

VRITELNC 'CNnR NUMBER OF STEFS IN 1UN6E-EUTTA NETROD'); 

lEADLN; 

lEADCN); 

VRITECEFFINDX.N); 

VRITELNC '); 

VRITELNC STEPS IN RUNGE-EUTTA*' ,H;S): 

ENOCMF*) 

ELSE 

C* READ FILE *) 

BEGIN 

RESET! EFFINDU: 
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iKCT(trriND2): 

miKcrriwi.NUKVAiis}: 

ICA0<Cr?lH01,OIVlCNCTM); 

roi i:>( TO Ntmvilis*! oo iuo(crriNoi.isPACcnn: 
roi I:>l TO HUnVAlLS DO IEAO(ErriMDI,TVni)i 
roi 1:>1 TO NURVALIS DO IEAD{EFriKDl,QVtn)i 
IEAD(EmHDI,HAXEFF)i 
ICAO(EFriNDl.Tdn]>; 

TVCI]:*T0tl]; 

ICAOtZFFINDl.POnil: 

1EA0(EFF1ND2,K); 

IND(*CISE*) 

00: (• INDATA •) 


nOCEOUIE OUTOATA; 
lECIN 

VIITEIMICFFODTD.'IAMFIOV OUTPUT' :H); 

WIITSLNCEFFOUTD,* M: 

VIITElN(EFFOUTD,‘ 

VHTEINIEFFOUTD,' 

niniKIEFFOUTD, 'VALUES OF INPDT':IO); 

VIITEINCEFFOVTD.' 

WIITEOKEFFOUTD,' '): 

niTELNCEFFOUTD/INLET STACHITIOK TENPEIATVIE (I)>‘:!D,TOtl]:lS:ie); 

VllTElN(EFr0UTD,' INLET STAGNATION PBESSVRE (PAU';S0,POni:lS:lO); 

VIITELNCEFFOUTD.'NUHSER OF lAOIATION SHlELOS>*:!)l,I(UKVALlS:lS)i 
niTELNIEFFOUTD,' 

VIITELNIEFFOUTD. 'TEMPERATURES OF RAOIATIATION SHIELDS' :SB) ; 

FOR I:>I TO NVMVALLS 00 VIITELNIEFFOUTD. 'SHIELD TEMPERATURE (E)':31. 1:2, TVm:IS:I<) 
, VRITELNIEFFOUTD.' 

V1ITELN<EFF0VTD,'HEAT INPUTS TO BiOIATION SHIELDS' ;Sii) ; 

FOR I;«i TO KUXVALIS 00 VEITELNIEFFOUTD, 'SHIELD HEAT (V)' :3I.I:2. ‘■':l.QVni:lS:ll); 
UIITEINIEFFOUTO,' 

VIITELN(EFFOVTD,’PIPE LENGTHS 8ETVEEN SHIELDS' ; St) ; 

FOR I:>1 TO NUHVALIS«1 00 VRITELNIEFFOUTO. 'SPACE LENGTH (M)':3I.I:I.'>':l,lSFiCErM:tS:l«) 
VRITELNIEFFOUTD,' 

VRITELNIEFFOUTO, 'S0LU7.I0N SECTION SIZE (M)>' ;S0,DIVLENCTH:1S:10); 

VRITELNIEFFOUTD, 'N12IMUN SHIELD EZCKANGER EFFECTIVENESS-' :SI,MAZEFF:1S:1I): 
VRimNIEFFOUTO,' '); 

VRITELNIEFFOUTO,' 

VRITELNIEFFOUTD,' 

VRITELNIEFFOUTD,' 

VRITELNIEFFOUTO,' 

VRITELNIEFFOUTD, 'VALUES OF OUTPUT’ 

VRITELNIEFFOUTD,' '); 

VRITELNIEFFOUTD,' 

VRITELNIEFFOUTD, 'MASS FLOV (EC/S)-' :5S,MD0T;2S;1B>; 

VRITELNIEFFOUTD, 'PIPE DIAMETER (M)-' :SO,OIi[LOOPCOUNTl:30:ll}; 

VRITELNIEFFOUTO,' 

VRITELNIEFFOUTD, 'TOTAL PIPE LENGTH (H)-':SI,LT0TAL:2I: 10); 

VRITELNIEFFOUTD, 'SHIELD LENGTHS' :<0) ; 

FOR I:«l TO KDKyALLS 00 VRimNIEFFOUTO, 'SHIELD LENGTH (M) ' :3C,I:2. '-' ;1,LTEHPCI1:1S: 11); 
VRITELNUFFOUTD,' '); 

VRITELNIEFFOUTO, 'TOTAL HEAT TRANSFER (V)*';SB,aVALLLNUMVALLS«l]:3l:ll); 

VRinLNIEFFOUTO, 'SHIELD HEAT TRANSFER'; 10); 

FOR I;-t TO HUNVALLS 00 VRITELNIEFFOUTO, 'SHIELD HEAT (V)':3I,I:2, '>':1,0VALLII]:1S:1I); 
VRITELNIEFFOUTD,' '); 

VRITELNIEFFOUTD, 'OVERALL EXCHANGER CFFECTIVENESS»'.SB,EFF[NUMVALLS«11:I0:10); 
VRITELNIEFFOUTD, 'SHIELD EFFECTIVENESS' ; 40 ); 

FOR list TO HUNVALLS DO VRITELNIEFFOUTD, 'SHIELD EFF':30,I:2,‘«' ;l.EFFm:lS:ll); 
VRimNIEFFOUTD, ' 

VRITELNIEFFOUTD, 'EXIT STAGNATION TEMERATURE II)-' ;5fl,TIiniUMDm;20;ll); 

VRITELNIEFFOUTD, 'EXIT STAGNATION PRESSURE IPA)-' :Sfl,PICNUHDIV]:2S:10); 

VRITELNIEFFOUTD.^ 

VRITELNIEFFOUTD,' '); 

VRimNI'FOR AN OUTPUT OF SECTION FLOW VARIABLES, INPUT 1. OTKERVISE I)ffUT I.'); 

REAOLN; 

lEAOIPXET); 

IF PIET « 1 THEM 
lECIN 

VRITEIEFFODTD, 'LENGTH':?, 'MACH NUMB ER * :H, 'STAC TEMP H ,' STAC PRES’tH); 

VRITEI EFFOUTD , ' TEMP ' : 1 0 , ' PRESSURE ' ; 1 fl , ' HUSSELT' : 1 0 , ' REYNOLD NUN ' ; i S ) ; 

VRITELNIEFFOUTD, 'SECTION HEAT' ; IS. 'TOTAL HEAT'tlS); 

VRimNIEFFOUTD,' '); 

FOR I:-l TO HUKDIV DO 
lEGIN 

VRITEIEFFOUTD.Zm:l;S,' '-.I.RUlilO:!,' ' ;4,Tg[I];lD:A); 

VRITEIEFFOUTD.' ' :3,P0III; IS:S, ' ' :2.TII):I;4, ' ';2,Pm:l:3); 

VRITEIEFFOUTD, • ' :2 .XUII) ;l:3, ' ' :5,REm:lfl;S,' ‘:S, 0111:10:7); 
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VimiNiEErOUTD. ' ’;S,QZ{I1:1I:7»; 

CKD(*rOI*> 

END:(*ir*» 

VIITElNrrOR A PIOT rilE EKni I, OTMEKVISE EHTEl IVi 
IEAniH;lEAO(Tll); 
ir TAI ) I THEN 
lECIN 

HllTElM(*rOI PLOT OF STACXATION TEHPERATVIE EXTEI 1 .'); 

VIITElHrrOR PLOT or STACXATION PRESSURE ENTER Z.'); 

HIITELNCrOR PLOT OF HEAT HANSFER ENTER 
READER; READ (TAD; 

VHTELNCPOUTD.NUHDIV); 

FOR !:■! TO NUXDIV DO 
RE6IN 

CASE TAR or 

1: WR1TELN(POUTD,TO(I]): 

2: VRmLN(POUTO,P«m>; 

3: VRlTELN(POVTD,aZ[Il)i 

END(<CASE«) 

END: (•FOR*) 

FOR !:■! TO NVXDir 00 VRITELNCPODTD.ZnJ): 

END;(*ir>) 

VRITELNCNUN8ER OF INTERATIONS>‘ ;30,LOOrCOUKT:S); 

VR1TELN(‘F0R INTERATIOH INFORHATIOK, INPUTl. OTHERVISE INPUT I* I: 

READLN; 

REAO(rXET); 

IF PREf ■ 1 THEN 

FOR I:>l TO LOOFCOUNT 00 

VRITELNI EFFOUTD . ’ INTERATIOlb ' : 21 , 1 ; 3 , 'DIAMETERS ' ; 21 .0 lAE 1 1 ; IS : II ) : 
END;(* OUTDATA «) 


FUNCTION FVR(N,P:RCAL):REAl; 
(* SOLVES FOR K>*P *) 

TAR A:REAl: 
lECIN 

A:sP*LK(N); 

PVR:«EZP(A): 

EXD:(* PVR <) 


FUNCTION R (N : REAL ): REAL; 

(t SOLVES OFTEN OCCURRINC FUNCTION D 
RESIN 

l:st.lf(E-t.l)(SQR(N)/2.l; 

EXD;(< I •) 


FUNCTION TEHP(TU,K:RCAL);fiEAL: 

(t CALCULATES LOCAL TEKFERATUHE *» 
lECIN 

nNP:>TO/S(K); 

END:(* TEMP I) 


FUNCTION PSCS(PO,N;ICAL):RCAL; 

(• CALCULATES LOCAL PRESSURE <> 

VAR A.C:REAL: 

lECIN 

C:>I(M); 

A:>ZP(1-R); 

riES:tPI*FVR(C,A); 

EN3;(i PRES *) 


FUNCTION DENSE(P,r;lEAL):REAL; 

(• CALCULATES LOCAL OENSITT <) 
lECIN 

DENSE: sP/tRtT); 

ENO;(* DENSE *) 


FUNCTION HU (T: REAL}: REAL; 

VAR A.S,C,0:1EAL: 

RECIN(*IIU»> 

A:>1.1041?2AS531E-D; 

l:sl.U«t77SM(<IE-l7*T; 

C:t.3.222>^t272m2E<l|tSOR(T): 
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c. 


0;>4.7IS700]frifm>l3«T*SQI(T}i 

OOi(*KV*) 


rmCTIOX COND(T;ICAl);RCAli 
ni l.I.C.OllEAL; 
IICIN(«CON]}UCTmn caicvlationm 

A:>I.I1S721II71; 

l:ilJ4KS2St.0MT; 

ff407S21M7*SQI<T): 

0:«4.S2H2l7<E-IO«T(SQI(T)i 

C0KD:>A*l«C«0: 

DO:(«COIID«) 


rUKTIOX rRANDTL(T;IEAL);tCALi 
VII A,l,C,0,E,F:REAL; 

IECIH(<PRAHDTl HUKBEI CALCULATION*) 

A:>I.S70232774; 

l:«I.IOSSI72!71*Ti 

C:a-7.SI324S32E-0S*SQR(TI; 

D:a4.32720S01C-«7*TiSQR(T): 

E:a-l.n«144SE-07*Sai(T)iSQI(T); 

r:at.l7403l42E-12<T>SQR(T)iSQI(T); 

PIAXDTL:>Atl4CtO«E»F: 

EXDi(iplANDTL*> 


FDXCTIOK )0UT(H1 ,N2,01 ,02 .REAL) tlEAL; 
(• USED IN CASE 1 OF NEVTOH *) 

VU A,C:1EAL; 

ICCIN 

A:-<Itl.f)f(2.0*(I-(.a)); 

C;>I(R2)/B(NU: 

HOVr:>(HI/H2>*FVB(C,A}~SOI(O2/0n; 
INDi(* }OUT «> 


FUNCTION 0NQUT(H1,H1;1EAL};BEAL; 

(* USED IN CASE 1 OF NEVTON *) 

VAI A.CilEAL; 
lECIN 

A:>I(M2)/B(M1): 

C:><X4l.O)/(2.0*(E-1.4}); 

OKOBT:««(l4l.8)/2.()*<B(NmB(H2))-1.0/SOI(N2))*HltPVRlA.C); 
EKD;(« OKOUT *) 


FUNCTION NFRIC(N1,K2,0,F,L:1EAL};REAL; 

VAB NFI:REALi 

. (* USED IN CASE 2 OF NEVTON >) 

BEGIN 

HFl:><SQR(K2)-Sai(Nl))7(EtSQR(Hl«N2))4(X«l.l)/(2.S*l)*LN(l(H2)(S0R(NI)/(l(Hl}«SaR(M2>))-4.l*F*L/O; 

KFIIC:<NFR; 

END:(* NFRIC *) 


FUNCTION DHrilCINl ,R2 :REAL) '.REAL; 

(* USED IN CASE 2 OF NEVTON •) 

VAI DNF:IEAL; 

BEGIN 

DHF:«2.l/<C*N2iSQI(N2))4(E4l)/E*H2*Sai(l(Hl>}/(PVR(Hl,4)<>PVB(l(N2).3)); 

OHFIIC:xOHF; 

IF DHT < 1.274 THEN (*FUNCTION UNSTABLE, KACH>1 APPROACHED*) 
OHFIIC:al.l: 

END;(* OKFRIC •) 


FUNCTION KrHST(N.O,HDOT,PI,TI :IEAL) :REAl; 

(t USED IN CASE 2 OF NEVTON •) 

VAR A,C:mL; 

BEGIN 

A:>B(H); 

C::(Xil.0W(2.1*(l.l~E)); 

NFIST:i(N*F0*PI>SQS(O)/4.S)*SQlT(I/(R*TI))*PVR(A,C)-)!!)OTi 
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00; (• liriST •> 


rmCTION DHriSKN. D, r I .TO :IEU) ; lElL ; 

(• osio IK cisc s or nevton *> 

VII l.c,t:k 2 lli 

ICCIN 

l;>l(H): 

!:■(]. 1*1-1. 9)/(2.IM1.I-E)): 

mnsT:>(ri*SQi(0)/i.i)*ro«sQiT(K/<o*Ton*(vvi(i.cMQi(H)«(K*i.i)/].o*rvi(i,E)i 

INO(* ONflST t)i 


rUKCTIOM MCVrOKdCT.’IKTECEl; Rl, 0 .r.l,LTOT,rO.TI;ICU):ICIti 
VU KKCV.XOLD,Dirr,OI,D 2 ,A,n}OT.OHr,Mr;IElLj 
eOOHT: IKTECEl: 
lECINt* NEVTOH *) 

Dirr:«i.ii 
cisi iiT or 
1 : lECIN 

(* Clicvims EXIT HICK NUHICI fOI SECTIOK WITH IIEI CHANCE *) 
OI:iO; 

02 :>r: 

iniO:>Sai(Dl/D 2 )*Nli 
WHILE (Oirr }> I.IOt) DO 
lECIN 

KKEV; «ROIO-HOUT(K1 .ROlO ,01.02)/ (ORODTtRl .ROLO) ) ; 
Oirr:*ABS(ROLO-KNCV); 

ROlOriRNEV; 

END(* WHILE *) 

EKD;(* I •) 

X: lECIN 

(< CiLCOLATES EXIT RACK rOE SECTION WITH ERICTIOK ONLY <) 
HOLO:il. 01 *Rli 
COUNT: > 1 : 

WHILE (Dirr )> f.OOOIS) ANO (COUNT ( II) DO 
lECIN 

ORE :«ONrilC(Rl .HOLD): 

Rr:>HriIC(Rl,H 0 E 0 , 0 ,r.l); 

RXZV:>R0l0-MrBIC(Rl,H0l0,0,F,L)/0Rr; 

Oirr:sABS(HOLO-HNEW); 

ir OHr > 1.0 THEN (*rUNCTION UNSTABLE, EXIT lOOP*) 
lECIH 

Dirr:>f. 10001 : 

HHEV:>t.l; 

CND:(*ir*) 

COUNT:>COIINT«t; 

HOLO:*HNEV; 

END(* WHILE *) 

END; ((CASE 2*) 

I: lECIN 

(< CALCULATES INLET RICH KUKBEI *) 

HOOT:iL: 

N01D:<R1; 

WHILE (oirr )« i.oioit) oo 

BECIN 

HNCV::KOLD-HrRST(HO L 0 , 0 , RDOT.ro, TO)/DRriST(KOLD.O,rO. TO); 
Oirr:>ABS(ROLD-KNEV); 

ROLO:sRNEV; 

END(* WHILE «) 

END(* 2 *) 

END:(* CASE t) 

HEVIQN:«ROLD; 

END:(« NEWTON «) 


PIOCEDUIE INLETKICH(D.KOOT,rO,TO:BEAL: VU R,T.r,BKO:BEAL) ; 
(* CALCULATES INLET CONDITIONS *) 

TAR A:IEAL: 

IECIK(* INLETHACI >) 

K:>NEVr0NO,H,0,0,HD0T,l.l,r0,T0); 

T:«nRP(TO,R); 

P:>riES(FI.R); 

IRO:xOENSE(P,T); 

END:(* INLETRACH >) 
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riOCIOVIC imCHAHCE(DI,Ofi'.TO,;i;REALi»l H.T.P.lHO.IE.r,NU.»:IEU); 
(< roi DUHETEI CHANCE *> 

HI A:IEAl; 

ICCIN 

N:aNCVr 0 M(l,H, 01 . 0 M,l,M>: 

T:>nHP(T 0 ,H); 

P;>PXE 5 (M.M): 

IIO:>OENSE(P,T>; 

IE:«IIi 01 /DXi 
IP IE (• 2101 THEM 

lECIN <« lANINAI riOV M 
P:>U.lflEi 

Eio(tirti 

ELSE 

lECIN (t TUIBULEKT FIOV *t 
P:al.Hf/PVX(IE, 0 . 2 ); 

IU:>l. 013 tpVR(IC, 0 .IMPVI(PI,M); 

EXD:(iELSEt) 

00;(* AIEACHAKGE •) 


PIOCEDOIE rXICTION(TD.DIA,F.L:IEAL;VAl K,T.P,IHO,PI.Q:IEAL); 
(• FOl FIOV VITH FRICTION ONLY •) 

VAI A.C,NEVH:tEAl; 

ICCIN 

XCVH: >NCVT0N(2 ,N. 0 1 A , F , E , LTOTU, I . I ) ; 

IF NCVH > t.l THEN NEVN:ai.O; 

A:>I(MCVNWI(N); 

C:><I«l.l)m.D<(I-I.O)); 

, PI:*PO<HINEVK*FVI(A,C): 

HrmNEVN; 

T:»TEHP(Tfl,N)i 

P:sPIES(PS,R); 

IHO:>OENSC(P,T): 

Q:<i 0 : 

EN0i{*FllCT10N*> 


FUNCTION DIAHETEI{HDOT,TO,Pe,K:SEAl);ICAL; 

(• CALCULATES OIANETCR CIVEN MOOT AND HACR •> 

VAI A,C:IEALi 

ICC1H(*0IAHETEI<) 

A:«4.R«NDOT*SQlT(R*TO/E>/(PIfN*PI)i 

DIAMETEX:a5QRT(A*PVR(B(H) C)I; 
ENDi(*DlAKETCR«) 


FIOCEOUIE SCCTIOX(HUKrOT:lKTECER; LENGTH, OIVL; REAL; VAI I:OIVISIONS; TAR NVNDIVrlNTECEl); 
(« DIVIDES PIPE LENGTHS INTO SKALLEl SECTIONS <) 

VAI 

ICSL,1EAL0IV:REAL: 

I.JiINTCCEI; 

BEGIN(*SECTtON«) 

IF NUHTOT « 1 THEN 
ICGIH 
KUHrOT;>0; 

END;(«IFt) 

(• CALCULATE RESIDUAL SOLUTION SECTION lENCTl t) 

ICALDIV:iL£KCTH/OIVL; 

NUNDIVisTRUNCIREALOIV); 

ICSL:«(ICALDIV-NVKDIV)*DIVL; 

IF ICSl ) I THEN 
ICCIN 

NUI(01V;<NVHDmi: 

END(<IF<) 

ELSE 

BEGIN 

ICSL;aOIVL; 

END;<*ELSE*) 

<• DEFINE lEKCTI SECTION ARRAY 
NVHDIV:>HUNTOTiNUNDIT; 

FOR Ii.NUKTOTd TO NUHDIT DO X[I]:.XII-13*0IVL; 

XIKVMD I VI ; 1 1 INUHO I V- n 4 RESl ; 

I£NUKDIV*IJ;alINUHDIVUOIVL; 

IKD;(<SECTIONt| 
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rRiC7:0K STlCTCKr(TI.TV,NV.Kr,RDOT.L:aill);IEAL: 

(• CUCULATCS MCV STACHATION TCHFEIATVSC tti riOV VI7H HEAT HAKSFEI O 
IICIN 

STACTEHP ; :TV» ( (Tl'tVM Cir ( (MUtF 1 1 KE * I ) / (fU07« CP ) ) ) ; 
nO:(*STAGTENP() 

rimCTlOK inmC(K,TI,TV,Plt.NV;IEAl: KEr:IXTEGEI);KEAl; 

(t SOLVES AK OFTEM OSED EUNCTIOH IN PICCSOffRE HEATTIANS *1 

VAI A,C,0,E:I£AI: 

lECIN 

C:>M.l4l«SQI(M>)/Tti 

IP IIT > 1 TKEX 
lECIN 

ua(iir*) 

ELSE 

lECIK 

E:iiU.A*PI/(KU>(1V-TI)}i 

END;(«ELSE*} 

irVNC:>AMD*ElC): 

END; (3 IFVKC •) 


PIOCESUIE REATniANS{TV,l.O,RDOT;lEAL; K:IKTECER; VAI H,TO,PO,T,P,lHO,RE,r,a,aE,KV;IEAL) 
(* CALCULATES ELOV TAIUILES FOl FLOW WITH FIICTIOM AKD HEAT TIAKSFSI <) 

VAI 

Cl.C2,Ca,C1,SnP,A,C:llCAL.; 

Tll,TA2,Kl.H2,FliEF,2]ELX,E;IEAL; 

lEMrlKTECEl; 

lECIN 

SEI:.I; 

(« FLIC EEI IF FLOV TUIBULEKT *) ' 

IP KU > 1. 14 THEN EET:>i; 

(< CALCULATE MEV STACKATIOX TEMFEIATURE <1 
PI:iPIANDTL<T); 

EP:tCOND(T); 

T01:*T0; 

TI2 :»STACTEHP(Tll,TV,HU,EP,MDOT,lJ ; 

(« STEPVISE SO urn OK USING RUNCE-KUTTE KCTHOD •) 

(• CALCULATE STEP SIZE *) 

STEP;sAISUT02-TOI)/N}i 

K1:«H; 

I:>1; 

WHILE I {> K 00 
lECIN 

CI:<RFUNC(N,TO.TW,FI,KU,KET>; 

IP Cl < 0 THEN Cl :*0; ((FUNCTION UNSTABLE*) 
C2:<IFUNC((N*CltSTEPr2.0),(T0>STEF/S.9),TV,PR.NU.XEr); 

IF C2 ( I THEN C2:>0;< (FUNCTION UNSTABLE*) 
C3:>IFUNC((NiC2*STCP/2.9),{T94STCF/1.9),TV,FR,NU.XEY); 

IF C3 < 9 THEN C3:s9;( (FUNCTION UNSTABLE() 
CA:»RFUHC(HtC3*STEF,TB*STEF,TV,Fi,NU,IEI); 

IF CA < 9 THEN Cl: *9; ((FUNCTION UHSTABl£() 
OCLN:>STEPf4.0((Clfl.|(C2«2.0(C3»C4); 

H;>K40ELN; 

TI:iTl4STEF: 

1:-Uli 

IF R ) I THEN ((FUNCTION UNSTABLE, EXIT LOOP*) 
lECIN 
I:>N«1; 

H:>l.li 

CND(*IP*) 

END;(* WHILE () 

IF (Cl m I) OR (C2 > I) OR (C3 < 9) OR (Cl « 9) THEN 
((FUNCTION UNSTABLE, EXIT LOOP*) 
lECIN 
I:aK4li 
R:*l; 

ENDi(((F*) 

H2:«N; 

A:>I(H2)/B(N1); 

C:s(K*1.9>ni.9((X-l.l)): 

PI:«P0*(N1IN2)(SQRT(T02/TBI}(FWR(A,C)! 
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T:«TCHP(T01,)»); 

uo:-ocMsc(r.T>: 

II:«4.l<KD0T/(PIt0*HIMT>); 

IF 1C <■ KOI THEM 

■ICIM (* LiHlNAl FLOW 

r;iU.I/IC; 

muMi; 

IND(*IF«» 

lUC 

IICIK (• TURIULEKT flOV M 
F:!il.m/PVIH1E,I.2); 
IU:il.ll3tpV]l(RE.0.l><rVI(PI.I.1)i 
CXD:(«EISE*) 
a:sCP*KDOTt(T02-Tin: 

QI:>QS*Q: 

m;i* lEArriAKS ■) 


PROCEDURE SCTROUKD(nr:IKTECEI: VAI Hit, lOVI.OIA: REAL) : 

(< 117ERS OIAHETEl SO THAT HACK tiLT U INTEIATED IK THE HAIM PROCIAM TO REACH OMC ■) 
ICCIM(<SETI0VH0«) 

IF EC! ■ 1 THEM 
ICCIN 
III:sOIAi 
EKD(*IF<I 

ELSE 

lECIM 

lOVI;>CiA; 

ENO;(«ELSE«> 

01i:>(l0VI«KII)n.l: 

EXOi(*SCTIOUNDt) 

PIOCESURC SVAP(VAI A, I : REAL): ' 

(t USED IT PROCEDURE BSORT 
TAR T:IEAL: 

IEC1M(*SVAP*I 

T:-A: 

A:-l; 

l:.T; 

EKDi(<SVAP*) 


PROCEDURE ISORTCSTART, TOP: INTEGER; VAR ARRT;COUNT; VAR CKT.-ICOUKT); 
(t USED IK FUNCTION HDOTHIN TO FIKO OPTIHUR MASS FlOV M 

VAI 

RCKT: COUNT: 

IKFECER: 

SWITCH :BOOlEAN; 

IECIN(t|SORT*> 

FOR I:iSTART TO TOP DO RCNTni:sCKTtn: 

REPEAT(< UNTIL SORTED VITH SMALLEST AT TOP M 
SVITCR.'<FALSE: 

FOR I: .START TO TOP-1 DO 
lECIN 

IF ARRim > ARirilUl THEN 
lECIK 

SVAP(AERTm,ARRTn*n); 

SVAPIRCKTIIl.lCXnun): 

SVITCH.v.TRUE; 

END(<IF*> 

EKD:(*FORt) 

UNTIL SVITCH ■ FALSE; 

FOR 1:«START TO TOP DO CNTm:>ROUND(RCKTIl]); 

END;(t|SORT>) 


FUNCTION KDOTMINIQ, TV: COUNT; T01,EFF:REAL; X: INTECER) :IEAt: 
(t OETERHIXCS OPTIMUM HASS FLOV O 

VAI 

QTOT.TOKAX.HDOT: COUNT; 

INDEX: ICOUNT; 

I. NCOKT: INTEGER; 
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IICIN<*H90TNIHM 
(• INITULIZI IKOCZ nUJ *) 

roi I:>l TO N 00 IKDEXm;*ti 

(• CUCVUTE SHIELD EXIT STACKITIOX TEI9EIJITVIES •> 

Tiiunij;«m; 

roi l:>l 70 M 00 TORAEm:«Err*(TVm-TOKiX[M})4TOIUEn.il; 
(« CIICUIATE HASS ELOVS AS DETEIMINEO IE SECTIONS «) 
OTOnil;>l; 
roi I:-l TO M 00 
lECIK 

OTOT[n;«QTOT[I>U4ami 

mm 1 1 : ■QTOTI 1 ] / ( CP* (TOKAXn 1-TI 1 n ! 

CW;(*roi»» 

(• SEilCH roi lARCEST HOOT «) 

ISOlTd.M, HOOT. INDEX); 

RDOTHIK:>n>OTIK]: 

CND:(«KDOTXIN*l 


I 

I 

r 


nOCEOinC IKREASSICN; 

(*0EF1NES EALUES fOl USE IN PAIAHCni LIST •) 
IECIN(*I)aCASSICN*) 

Pli:>Pim; 

TiraTEll: 

Pl:>pm: 

UOl.-iXROCIl; 

Pl:»fn]; 

. IEl:>lcm; 

01:a0m: 

axE;>oxm; 

EXD:<*INlEASSICKt) 




i 

I 


nOCEDOIE OUTHEASSICN; 

<> DEFINES NEV VALUES OBTAINEO FXOH FAIAHEni LIST >) 
IECIN(*OUTlEASSIG)i*> 

I:*Ut: 

Nm:«Nl: 

TSm.-TOl: 

FIIII:>F9I; 

Tni:.TI; 

Pm:>Pl; 

IH0m;«IH01; 

FEUkFI; 

ICII]:«iEl; 

Qm:«Ql; 

axm:*QXl; 

NVm;«NUl; 

ENOi(tOUTlEASSICNt) 


t riOCEDURE INITIALIZE; 

(• FREFASES VARIABLES FOR USE IN HAIN PROCRUf *) 

. BECIN(» INITIALIZE*) 

(‘OETERNINE HINIHUH MASS FLOW*) 

, HDOT:»KIJOT«IN(QV,TV,T»(n,MAXErF,HUHViLlS); 

(• DEFINE riBST SECTION SIZE •) 

. XCll;sOIVLCNGTX; 

(• DIFIXE IMACINAHT LAST PIPE LENGTH «) 
LSPACC[NimVALLS42];>D; 

1 («STABT LOOP COUNTER*) 

LOOPCOOHT:il; 

((INITIAL HACK NUMBER GUESS*) 

Nin:<ij: 

. (* INITIAL DIAMETER GUESS. ASSUMES INLET XACH OF «.l *) 

DIAm:>OIAHETEB(MDOT, Toil), FOCII. Mill): 

. <• CALCULATE KIHIMUM DIAMETER, )'ACH « 1 *) 

lOVBOUND;*DIAMETER{KDOTJ9m,Pflm,I.«); 

(* DEFINE XAIIKUM DIAMETER *) 

HIBOUNS;<tOO*lOVBOVKD; 

END; ((INITIALIZE*) 


PROCEDURE REINITIAL; 

(* PREPARES 7ARIABLES FOR USE IN MAIN PROGRAM *) 
BECIN(*RCINITIAl*) 
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(•INlTIAUZt COUKTIISM 
I:at< 

rirccirrs.is 

VKLCirr:>t; 

SMCECKTiit! 

(iiNCicHiHr loor covnth*) 

ioorcouKr;>ioorcouKT«ii 

(• OiriNC HIST tin tCNCTI () 

iTom:*LsrACEnii 

(tCILCUllTE ENHANCE COKDITtOK!*) 

■En]:«4.0*NOOTf(PI*l10(Tltn)*DlA(lOOFCOUKr]); 
ir 11(11 <■ 2018 THEN 

lEClN (• LIHINAR ElOV M 
rniuu vimii 

END(*ir*) 

EISE 

IE6IN (( TURIUIENT ElOV M 
fn]:sO.I4l/FVI(RE[|],3.2); 
aiUCn:>l.l23<FVR(REm,O.I)«FVR(FRANDn(TOCn}.l.4):. 
ENO:(«ElSEi) 

(tCAECVUTE INLET CONDITIONS*) 

INIETHACHIDI At LOOFCOUim . MOOT, FOn 2JIC 1 1 ,Hm ,TI 11 , FC 11 ,IH0C in ; 
(t INITIALIZE COKHOL CARD*) 

HnST:>l; 

END;(*IE1NITIAI*) 


lECIHI* MAIN FHOCRAH *) 

. VRITEINCTO ENIER DATA, ENTER I. OTKERVISE, ENTER I.'); 

REAOLN: 
lEAO(TAI) ; 

INDATA(TAI); 

INmALIZEi 

REPEAT (* DNTIl THE FLOV IS CHOEED *) 

REINITlALi 

WILE (SPACECNT (« NUNVALLS«1) AND (HTEST < 2) DO 
lECIN 

(* TEST FOR CHANCE IN TENFERATURE BOUNDARY *) 

IF Zm ( (LTOTAL-OIVLENCTHIZ.n THEN 
EECINI* FRICTION FLOW *) 

(* OIYIOE FRICTION LENGTH INTO SECTIONS *) 

SECTION! I, LSPACEl SPACECNT] .DIVLENGTH, E.NUKDIV) ; 

<* DEFINE END HACK NUHBER FOB TEST AT END •) 

HCNUHDIV1;«I; 

VHiLE a ( maw and (ktest < ?> do 

BEGIN 

L:>ZllM(M]; 

INREASSICN; 

FRICTtONITOl ,DIAtlOOPCOUNT] .FI ,L,H1 .T1 ,Fi ,AHOI ,P01 .01) ; 

OUTRE ASSIGN; 

IF Nil 1 ) 1.17 THEN KTEST: >2 (*EZ1T LOOP*) 

END:(*VHIIE*) 

SPACECNT; >SFACECNT*1; 

FIFECNT.-xPIPECNTti; 

END!* FRICTION FLOV *) 

ELSE 

BEGIN!* 7EKFERATURE BOUNDARY CHANCE *) 

QVALLIVAlLCNT];xO; 

LTCKFCVALLCNT]:sl; 

EFFCVAllCKT]:>l; 

T0SAVE:>T0m; 

VHILE (QVIVALlCNn ) OVALLIVALLCNTn AND lEFFIVALLCNTI ! HAZEFF) AND IHTEST < 2) DO 
BEGIN 

INREASSICN; 

HEATTRANSITVIVALICNTI .OIVLCNCTK.OIAILOOPCOVNT] .NDOT.N.Hl .TOl.FBi.Tl.Pl .IHOl .REl . 
OUTREASSICN; 

QVAllIVALLCNT] : lOVALL tVALLCNTI «Qi ; 

EFFIVALLCNTl : X ( TO t n -TOSAVE ) / ! TV! VA LlCHTl-TO SAVE ) ; 

LTEKF CVAILCKT] : xLTEKP CVALLOTT] f 0 1 VLENCTH; 
imrxZCMnOIVLENCTH; 

IF Hill ) 0.17 THEN JrrEST;x3;!*EZIT LOOPS*) 

END;!tVHIlE«) 

LTOTAL:>LTOTALfLTEKP(VALLCNr]FlSPACElSPAC£CNTl; 

FIPECKT:rpiPEC)tT*I; 

VAllCKT:xVALLCHT4l: 

END;!* TEKPERATVRE BOUNDARY CHANCE *) 

END:!*VHILE<; 
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ORIGINAL PAGE 13 
OF POOR QUALITY 


DUtifOiiiioopcouimi 

SCnOtfKO ( RT eST. H I lOURO . to VBOVMO , 8 un ; 

oucioorcoimT»ii,*«D(M; 

8 imL (MtNVKDIVl > I.f 7 ) AND (PlfCCXT >• 2 *NUf{VAllS»l) 01 (LOOPCOUNT } 31 ); 


(tCUCUUTI OVUklL EmCTIVCNISS *) 

irr(NVNVAiis« 1 1 :>(T0 ()(ura)iv]-Ta( 1 1 ) / nvonnwAiisi-Tin i ) ; 

(• CUCULATE TOTAL HEAT TAANSEEI «) 

WALlCKU)(VALtS« 1 ] : «m>0T^CP* (TO(NimOIVl-TO( 1 ] ) ; 

VimiNt 'OUTPUT STOIED IN FILE EFFOirrO.'); 

VIITEINC INPUT STOIEO IN FILES EFFINDl AND EmNDl.'); 
OUtDATA: 

EAD.C1NA1N PIOMAH*) 


ORIGINAL PAGE 
OF POOR QUALITY 


eiNO. 

TO CNTCR DATA, ENTER i. OTHCRVI SE , ENTER 0. 

? I 

INPUT DATA AS PROMPTED BY TERMINAL 

WHEN ASKED FOR A SERIES OF DATA SUCK AS THE SHIELD 
TEMPERATURES, INPUT THE DATA Hi THE ORDER FROM DEVAR 
TO THE OUTER BOUNDARY. 

ENTER THE NUMBER OF RADIATION SHIELDS 
? 3 

NUMBER OF SHIELDS* 3 
ENTER SINGLE DIVISION LENGTH (M) 

r 0.002 

DIVISION LENGTH (H)- 0.00200 

ENTER PIPE LENGTHS NOT ON SHIELDS (H> 

1 O.Olf 
t 0.021 
1 0.024 
t 0.02A 
T 0.034 
? 0.035 

ENTER TEMPERATURES OF RAO I AT 1 01^ SHIELDS (K) ■ 
t 21.7 
? 3f.7 
? 45.4 
? 101.1 
? 144.4 

ENTER HEAT INPUTS OF RADIATION SHIELDS (V) 

? 0.08B4 
? 0.054 
f 0.0714 
?> 0.1384 
? 0.8242 

ENTER MAXIMUM HEAT EXCHANGER EFFECTIVENESS AS A FRACTION 
*0.98 

SHIELD EFFECTIVENESS* 0.7800000 

ENTER INITIAL STAGNATION T (X) 

? 4.23 

INITIAL STAGNATION TEMPERATURE (X)* 4.23000 

ENTER INITIAL STAGNATION P <PA> 

? 2100.0 

INITIAL STAGNATION PRESSURE <PA). 2100.00000 

ENTER NUMRER OF STEPS IN RUNGE-KUTTA METHOD 

? 4 

STEPS IN RUNGE-KUTTA* 4 
OUTPUT STORED IN FILE EFFOUTD. 

INPUT STORED IN FILES EFFINDl AND EFFIND2. 

FOR AN OUTPUT OF SECTION FLOW VARIABLES, INPUT 1. OTHERVISE INPUT 0. 
t 1 

FOR A PLOT FILE ENTER 1, OTHERWISE ENTER 0 
? 0 

NUMBER OF INTERATIONS* 13 
FOR INTERATION INFORMATION, INPUTl . OTHERVISE INPUT 0 
f 1 

2.804 CP SECS, 102714B CM USED. 

/ 
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niuis or iKPVT 


iniT SnCKiTION TtKPtlATUIC (K)> 4.2300I0I9H 

INICT STAGNATION PIESSUIG (PA)<il 10 .1810010101 
NUHICI or lADIATION SHICLDSii S 

nHPEIATUItS or RADIATIATION SHIELDS 
SKICID TENPEIATUIE (I) 1« 11.7800300000 
SHIELD TENPEIATUIE (1) 2> 37.7000101000 
SHIELD TENPEIATUIE (E) 3> 03.4008000000 
SHIELD TENPEIATUIE (E) 4> 101.1100001011 
SHIELD TENPEIATUIE (E) S> Mi.OIOIOIOOlO 

HEAT INPUTS TO lADIATIOH SHIELDS 
SHIELD HEAT (V) 1> 0.0104000000 

SHIELD BEAT (V) 2> 0.0340000011 

SHIELD HEAT (V) 3* 0.17UOOOOOO 
SHIELD HEAT (V) 4- 0.1300000000 

SHIELD HEAT (V) 3> 1.1242000800 

pip: LENGTHS lETVEEN SHIELDS 
SPACE LENGTH (H) 1> 1.0170000000 

SPACE LENGTH (N) 2> 0.0210000000 

SPACE LENGTH (H) !■ 1.1240000010 

SPACE LENGTH (N) !■ 0.0200000010 

SPACE LENGTH (H) U 0.1340010001 
SPACE LENGTH (H) 0> 0.0330001010 

SOLUTION SECTION SI2E (N>- 0.0020103000 

HAIIHUH SHIELD EXCHANCEl EFFECTIVENESS- 0.7110800001 


VALUES or OUTPUT 


HASS FLOW (EG/S)> 
PIPE OIANCTEl (Hl- 

TOTAL PIPE LENGTH (N)- 
SHIELD LENGTHS 


SHIELD LENGTH (N) 1- 
SHIELD LENGTH (N) 2- 
SHIELD LENGTH (H) 3- 
SHIELD LENGTH (N) 4- 
SHIELD LENGTH (H> 3- 


0.0320001010 
0 1100080900 
0.088 0101000 
0.0040000000 
0.0310100000 


TOTAL HEAT TRANSFEI (V)> 
SHIELD HEAT TlANSFEl 
SHIELD HEAT (V) 1- 0.0704774301 

SHIELD HEAT (V) 2- 0.0442481133 

SHIELD HEAT (V) 3- 1.1081734341 

SHIELD HEAT (V) 4- 0.1377374111 

SHIELD HEAT (V) 3- 1.7714723447 

OniALL EXCHANCEl EFFECTIVENESS- 
SHIELD EFFECTIVENESS 
SHIELD EFF U 0.4133727431 


SHIELD EFF 2- 
SHIELD EFF 3- 
SHIELD EFF 4- 
SHIELD EFF 3- 


0.3047710717 

0.2703141314 

1,2834743320 

0.7143441371 


EXIT STAGNATION TENEIATUIE (X)- 
EXIT STAGNATION FIESSUIE (PA)- 


0. 9100014272 
1.0018244434 

0.3338001000 


1.1723328433 


0.707814C440 


143.1470237333 

447.2377023794 


LENGTH 

NACH NUNHEI 

STAG TEHP 

STAG PIES 

TEHP 

FIESSUIE 

HUSSELT 

lETHOLD NUH 

SECTION HEAT 

TOTAL HEAT 

0.10200 

0.02143443 

4.230800 

2100.80001 

4.2273 

2077.173 

3.441 

713.43327 

1.0810001 

1.8108080 

0.10400 

0.02147770 

4.230000 

2073.77130 

4.2273 

2073.143 

3.440 

713.43327 

0.8009000 

1.0080810 

0.10400 

0.02133723 

4.231001 

2071.73004 

4.2273 

2071.140 

3.440 

713.43327 

1.0000000 

1.0000008 


o”rooi QUALITY 



Mim 

Mlllllll 

1.230011 

2017,73127 

1.2273 

2007.123 

3. Ill 

713.13327 

I.IOOOIOI 

MOIOOII 

• 

Mint 

Mlllllll 

1.230001 

2003.72771 

4.2273 

2103,117 

3. Ill 

713.13327 

0.8000001 9c 

0.1000306 

I.OIOOOOl 


MUSI 

M11II377 

1.230001 

2077,73101 

1.2273 

2077.111 

3.110 

713.13327 

MOIOOII 


I.IHH 

1.11171117 

1.231001 

2073,73737 

1.2273 

2073.123 

3. Ill 

713.13327 

I.OOOOOOl 

1 OlOIOSi 


I.IUtI 

M1171711 

1.23IIII 

2071.73317 

1.2173 

2171.137 

l.lll 

713.13327 

I.IIIOOOO 

MOIOOII 


Mini 

M11717II 

1.23I0II 

2017.77077 

1.2273 

2117.137 

3.141 

713.13327 

MIOIOOI 

MIOIOOI 


Mini 

1.11113111 

1.23IIII 

2011.10711 

1.2273 

2113.111 

3.441 

713.43327 

I.IOOOIOI 

MOIOOlO 


l.lllll 

M137331I 

1.II17II 

2013.72107 

3. nil 

2142.711 

3.441 

472.71171 

1. 1172071 

1.0172071 


M21H 

Mlllllll 

1.711777 

2013.12137 

3.7137 

2042.273 

3.441 

137.14077 

MI7017I 

1. 1112347 


MI3II 

I.I171171I 

1.727117 

2113.17307 

1.7237 

2141.714 

3.441 

114.11137 

1.1141742 

1.0211331 

i ‘ 

Mini 

1.11712711 

7.311071 

2112.71072 

7.3137 

2011.211 

3.441 

S77.15<77 

MI44I77 

M27I211 


Mini 

M301III1 

1.211711 

2012.31470 

1.2122 

2141.733 

3.441 

333. 13771 

1. 0041171 

1.0313012 

t 

Mini 

1.13173377 

7.I2I1I1 

2011.71111 

7.1231 

2041.211 

3.141 

333.43237 

1. 1042737 

M103I37 

, 

Mini 

I.I331I11I 

7.711131 

2041.33771 

7.7311 

2137.141 

3.441 

311.54301 

I.IIIOSSI 

1.0141373 

4 

Mini 

M3131I7I 

11.131012 

2011.07114 

11.1217 

2037.042 

3.441 

177.31121 

M03I2II 

M3214I3 


Mini 

1.11112701 

11.172311 

2011.11031 

11.1177 

2131.141 

3.141 

112.23717 

I.IISS772 

I.03IIISS 


Mini 

I.I3I137II 

11.721341 

2010.12177 

11.7211 

2037.011 

3.411 

111.31224 

MI33I24 

1.1431211 


i.niii 

1.13737113 

12.332217 

2037.10171 

12.3213 

2037.211 

3.141 

134.12412 

1. 1031244 

MII3S17 


Mini 

I.I3I13I01 

12.711211 

2037.03171 

12.7037 

2034.313 

3.441 

111.73141 

MI1I7II 

1.1731133 

* 

Mini 

1.03707131 

13.111373 

2030.17311 

13.1337 

2133.071 

3.141 

131.13751 

1.0014343 

e. 0711028 


Mini 

1. 11711111 

13.713371 

2037.70373 

13.7711 

2053.111 

3.141 

123.17337 

1.1011231 

MI2S270 

' 

Mini 

1.01033121 

11.177171 

2057.27701 

11.1717 

2031.171 

3.141 

117.33103 

1.0111774 

1.0147214 


l.lllll 

1.01122077 

11.717111 

2031.117.(7 

11.7111 

2033.733 

3.441 

117.77271 

1.8837731 

0.1701 771 

• 

Mini 

1.01127111 

11.717111 

2033.77011 

11.7107 

2131.171 

3.141 

117.77271 

I.OOOOOOl 

1.0701774 


Mini 

I.I1I32I7I 

11.717111 

2031.31313 

11.7117 

2010.373 

3.110 

117.77270 

I.IOOOIOI 

I.I70I77I 

i 

1. 11711 

1. 01131213 

11.717111 

2011.11277 

11.7107 

2HS.714 

3. Ill 

117.77271 

I.OOOOOOl 

1. 0701771 


Mini 

I.I1113I7I 

11.717111 

2015.77211 

11.7107 

2013.012 

3.110 

107. 77271 

MOOOOOl 

1.1701771 


Mini 

1.01117117 

14.717111 

2013.30133 

11.7101 

2011.370 

3.141 

117.77271 

I.OIOOOOl 

I.07GI77I 


3. mil 

1.14131331 

11.717114 

2010.13133 

11.7111 

2037.701 

3.141 

117.77271 

I.IOIOOOl 

M70I77I 

§ 

i.mii 

1.11117771 

11.717414 

2037.77301 

11.7111 

2033.033 

3.141 

117.77271 

I.OOflIOII 

1. 0701771 

M- 

1. 14711 

I.I11I114I 

11.717114 

2033.31377 

11.7101 

2032.311 

3. Ill 

117.77270 

MOIOOII 

M71I77I 


1.14711 

1. 11170701 

11.717414 

2032.13101 

11.7111 

2027.713 

3. Ill 

117.77271 

I.IOOOIOI 

0.0701771 

f- 

M71II 

1. 11171377 

11.717114 

2027.77011 

11.7117 

2027.013 

3.111 

117.77271 

MOIOOII 

1.1701774 

J 

M72II 

I.011I1III 

11.717111 

2027.31347 

11.7117 

2121.314 

3. Ill 

117.77271 

i.nioool 

1. 0701771 

J 

M71II 

I.I41II3I3 

11.377714 

2021.12771 

11.3172 

2023.114 

3.111 

313.71114 

1.1137771 

1.1011717 


M74II 

1. 11111323 

11.137721 

2023.12010 

10.1111 

2021.032 

3.141 

343.01720 

1. 1133144 

I.117I1S3 

y 

M7III 

1.01107111 

17.111772 

2021.31371 

17.1111 

2020.113 

3. Ill 

317.32417 

1. 1121772 

1.1307123 


l.lllll 

1. 11711203 

21.113132 

2023.17312 

21.1271 

loiii.ni 

3.141 

331.73443 

1.0123101 

1.1131224 


l.lllll 

1.11132121 

22.312717 

2021.71730 

22.3227 

201M13 

3.441 

31I.IIQ3I 

1. 0111231 

1.1317177 

' 

l.lllll 

1. 01110017 

22.312717 

2010.03771 

22.3221 

2011.337 

3.141 

311.40151 

I.IOOOIOI 

1.1317177 


I.II4II 

1.03117711 

22.312707 

2013.71013 

22.3221 

2011.272 

3.140 

310. 11038 

O.IOOOOOl 

1.1317177 


l.lllll 

1. 03173131 

22.312707 

2012.10323 

22.3223 

2000.171 

3. Ill 

311.10031 

I.IOIOOOl 

1.1317177 


i.inii 

I.I11I37I2 

22.312717 

2007.11177 

22.3221 

2003.111 

3.141 

31I.II03I 

e.ioooool 

1.1517177 

1 

i.mii 

0.03171711 

22.312717 

2001.31023 

22.3221 

2002.131 

3.141 

31I.III5I 

I.IOIOOOl 

1.1317177 


.1.17101 

1.03177727 

22.312707 

2013.17017 

22.3223 

1771.733 

3.141 

31I.II03I 

o.ioooni 

1.1517177 

1 

1.17101 

1.03207732 

22.312707 

2000.10302 

22.3222 

1773.011 

3.141 

311.41830 

I.IOOOIOI 

1.1517177 


M7I0I 

I.IS21S711 

22.312707 

1777.33730 

22.3222 

1772.007 

3.141 

311.10031 

I.OOOOOOl 

1.1517177 

1 

1.11011 

0.03223777 

22.312707 

1771.27317 

22.3221 

1707.737 

3.141 

311.10030 

I.OOOOOOl 

0.1517177 


1. 11201 

I.IS231SS7 

22.312707 

1771.21111 

22.3221 

1701.141 

3.141 

311. 11031 

0.1000081 

1.1317177 


1.10111 

1.03237711 

22.312707 

1701.13111 

22.3220 

1703.101 

3.141 

311.10030 

I.OOOOOOl 

1.1317177 

] 

i.inoi 

0.15211013 

22.312707 

1703.07333 

22.3217 

1700.331 

3.141 

31I.IIC3I 

I.IOIOOOl 

1.1317177 

1 

1. 10110 

0.03I277II 

23.111313 

1703.27337 

23.1310 

1771.111 

3.141 

271.11117 

1.0210101 

1. 1130271 

.1. 

1. 11001 

0.11777173 

27.117721 

1701.20313 

27.0031 

1773.302 

3.141 

211. 03104 

1. 1273712 

1. 1103717 


l.lllll 

1.01303372 

32.233003 

1770.07031 

32.2121 

1772.311 

3.141 

230.31721 

1. 0217101 

0.2373773 


1.11100 

I.0II0I177 

33.327211 

1771.21272 

35.1757 

1717.177 

3.141 

233.11712 

0.I2S711I 

1.2431203 

I 

l.lllll 

1. 01117720 

33.327211 

1772.17017 

33.2733 

1741.772 

3.141 

233.11712 

O.IOOOOOl 

1.2131203 

1 

l.lllll 

I.III33I77 

33.327211 

1710.11012 

33.2732 

1710.703 

3.141 

233.11712 

I.OOOOOOl 

1. 1131213 


l.lllll 

1. 01117330 

33.327211 

1711.01710 

33.2730 

1751.017 

3.141 

235.11712 

1.0108001 

0.2131203 


l.lllll 

I.IIII112I 

33.327211 

1757.77372 

33.2710 

1732.732 

3.440 

233.11712 

0.0080081 

1.1131213 


l.lllll 

0.0II7S3I7 

33.327211 

1755.70307 

33.2711 

1710.111 

3.141 

235.11712 

1.0010080 

0.2431203 


l.lllll 

I.0III7372 

33.327211 

1731.03131 

33.2711 

1711.341 

3.140 

233.11712 

0.1008010 

1.2131283 

1 

l.lllll 

0.01703133 

33,327211 

1717.71121 

33.2711 

1710.173 

3.141 

233.11712 

1.0000000 

1.2131203 

* 

1.13001 

1.01717337 

33.327211 

1713.17123 

35.2737 

1731.371 

3.140 

233.11712 

1.0000000 

1.1131203 


1.13101 

1.01731737 

33.327211 

1737.12110 

33.2737 

1732.303 

3. Ill 

233.11712 

1.0800001 

1.2131203 


1.13101 

I.II71S777 

33.327211 

1733.35171 

33.2733 

1720.221 

3. Ill 

233.11712 

O.IOOOOOl 

1. 1431213 

t 

I.13II0 

0.01710277 

33.327211 

1731.10313 

33.2732 

1721.133 

3. Ill 

235.11712 

I.OOIOOOl 

1.2131283 

* 

1.13101 

I.II771I3I 

33.337211 

1727.11131 

33.2730 

1720.151 

3.110 

233.11712 

0.1000000 

1.2131203 


1.11011 

I.0I7I7OS7 

33.327211 

1723.31711 

33.2721 

1713.713 

3.141 

233.11712 

I.OOOOOOl 

1.2131203 


1.11101 

1. 17371137 

11.130312 

1717.31017 

11.3023 

1710.021 

3.140 

217.77271 

0.1335704 

1. 3111707 


1.11111 

0.07717717 

17.711271 

1713.02011 

17.1131 

1701.732 

3.141 

170.0112? 

1.0327311 

1. 3171217 

4 

1. 11101 

I.II1I3101 

33.772211 

1707.71000 

33.0130 

1170.331 

3. Ill 

171.17140 

0.IS12S7I 

1.1210777 


l.lllll 

0.OI170SII 

33.772211 

1701.20121 

33.0122 

1072.771 

3.141 

171.17111 

I.IOIOOOl 

1.1101777 

* 

l.llOOl 

I.IIS1S303 

33.772211 

1070.10010 

33.0113 

1007.232 

3.141 

171.17110 

O.IOOOOOl 

M20I777 


l.lllll 

1.01510137 

33.772211 

1073.11171 

S3. ns? 

1801.113 

3.141 

171.17110 

I.OOOOOOl 

1.1101777 

. 

1.11101 

0.0ISI37I3 

33.772211 

1007.10377 

33.0377 

1071.000 

3.141 

171.17110 

0.1008000 

1.1201777 


l.lllll 

O.OIS711II 

33.772211 

1002.03117 

33.1371 

1070.301 

3.141 

171.17110 

I.00I83SI 

1.1101777 


I.IIIOI 

I.III17207 

33.772211 

1071.10171 

33.0303 

1011.701 

3.140 

171.17111 

I.IOIOOOl 

1.1281777 


l.llOOl 

0.01113127 

33.772211 

1070.71101 

33.0373 

1037.277 

3.141 

171.17110 

I.OOIOOOl 

1.1201777 


I.IIIOI 

I.0III72S1 

33.772211 

1013.33101 

33.1317 

1033.101 

3. Ill 

171.17111 

1.1010008 

1.1201777 


l.lllll 

I.0II7S3II 

53.772211 

1037.71277 

33.0337 

1010.033 

3. Ill 

171.17141 

i.onoiii 

1.1201777 


ORIGINAL PAGE IS 
OF POOR QUALITY 



0.00722120 

53.772244 

1131.11041 

33.1331 

1012.115 

1.440 

171.47141 

0.1000001 

1.1201777 

1. 14100 

0.00714001 

53.772244 

1111.52135 

53.1312 

1034.744 

1.441 

171.47141 

0.0000001 

0.1101777 

0.17000 

0.0077SOSO 

33,773244 

1112.70372 

33.1331 

1131.101 

1.440 

171.47141 

O.OOOOOOlo^ 

0.1000000^° 

0.1201777 

o.mio 

0.00003037 

53.771244 

1137.27244 

33,1325 

1125.131 

3.440 

171.47140 

0.1201777 

0.I7000 

0.00030111 

31.7/S244 

1131.42701 

33.1314 

1117.757 

1.440 

171.47141 

0.0000000 

1.1201777 

I.1740e 

0.00030077 

33.772244 

1123.77117 

33.1301 

1111.044 

3.441 

171.47141 

1.0000000 

1.1201777 

0.17100 

0.00003734 

33.772244 

1120.30772 

33.0277 

1101.341 

3.440 

171.47141 

0.1000100 

0.1201777 

1.10000 

0.00711024 

33.772244 

1111.43111 

53.1270 

1102.431 

3.441 

171.47141 

1.0000000 

1.1200777 

o.iim 

0,07030000 

45.277411 

1004.47701 

45.0471 

1772.122 

3.440 

132.12477 

0.0754737 

0.5143534 

0.10100 

0.10734737 

74.101511 

1774.11721 

74.1114 

1777.447 

3.441 

134.77311 

0.1711437 

0.4107171 

0.10400 

0.11314331 

04.772432 

1701.74074 

14.4037 

1743.233 

3.440 

125.37113 

0.1175314 

0.7002721 

0.10000 

0.12204071 

74.737112 

1771.01115 

74.2321 

1710.717 

3.440 

114.77750 

1.0121440 

0.7027317 

o.moo 

0.12730520 

113.131112 

1733.12113 

101.0412 

1711.735 

3.441 

110.17722 

1.1717371 

l.lS417S:i 

o.ifloo 

0.13547032 

113.014347 

1737.31321 

112.3334 

1711.117 

3.440 

105.31717 

1.0412111 

0.7207177 

0.10100 

0..11I2S11S 

117.311122 

1717.71141 

111.7204 

1417.131 

3.440 

111.17517 

0.1517233 

1.7751130 

0.10400 

0.11435571 

121.717172 

1474.34351 

121.1277 

1444.371 

3.440 

77.11714 

1.1137371 

1.0211121 

0.10000 

0.15107121 

127.347417 

1473.72117 

121.3131 

1412.131 

3.441 

77.11177 

0.1374331 

1. 1301151 

0.20000 

0.13SSS411 

132.710121 

1417.71147 

131.7111 

1417.071 

3.441 

73.31711 

1.03!5(11 

1.0073112 

0.20200 

0.13702530 

13S.M1441 

1421.74151 

131.7317 

1571.421 

3.440 

71.15215 

1.12..;731 

1.1137374 

0.20100 

0.14377044 

131.200031 

1571.14107 

134.7443 

1341.134 

3.441 

73.21711 

1.1175707 

1.1335111 

0.20400 

0.14000373 

110.033137 

1371.13707 

131.7207 

1331.474 

3.441 

72.11171 

1.0153142 

1.1170414 

0.20000 

OJ72210S4 

HI. 171722 

1512.12114 

110.0170 

1505.277 

3.440 

71.71117 

1.1122235 

1.1412701 

0.21000 

0.17411017 

112.412551 

1513.31117 

111.1111 

1171.731 

3.440 

71.51414 

0.0075734 

1.1701137 

0.21200 

0.10071104 

113.177172 

1113.11340 

111.7141 

1113.711 

3.140 

71.20111 

0.0075015 

1.1743112 

0.21100 

0.10325471 

111.171317 

1132.32720 

112.3524 

1111.331 

3.440 

70,74571 

0.1031340 

1.1112112 

0.21400 

0.17000771 

111.730417 

1121.31171 

113.1011 

1371.115 

3.440 

70.77270 

0.0015407 

l.lllOOSI 

0.21000 

0.17305711 

115.117124 

1317.31200 

113.3231 

1311.227 

1.441 

70.44711 

0.0035170 

1.1723521 

0.22000 

0.20021327 

115.117124 

1333.12173 

113.2257 

1300. 7S1 

3.440 

70.44711 

0.1610001 

1.1723521 

0.22200 

0.20501370 

115.117124 

1311.15014 

113.1174 

1272.421 

1.441 

70.44711 

0.1010100 

1.1723321 

0.22100 

0.21207724 

115.117124 

1201.03411 

112.7753 

1231.733 

3.440 

70.44711 

0.0010101 

1.1723521 

0.22400 

0.21100023 

115.117124 

1211.31147 

112.1347 

1173.727 

1.440 

70.44711 

1.0000001 

1.1723521 

0.22000 

0.22412735 

HS. 117124 

1205.54157 

112.4717 

1155.141 

3.440 

70.44711 

0.0000001 

1.1723521 

0.23000 

0.23101111 

115.117124 

1145.32111 

112.5147 

1113.313 

3.440 

70.44711 

0.1000001 

1.1723521 

0.23200 

0.21132301 

115.117124 

1123.13311 

112.3011 

1047.331 

3.440 

70.44711 

1.0000000 

1.1723521 

0.23100 

1.25511371 

113.117124 

1077.72353 

112.0324 

1023.207 

3.440 

70.44711 

o.ieioool 

1.1723521 

0.23400 

0.24754717 

113.117124 

1033.13111 

111 .7501 

771.314 

3.440 

70.44711 

1.0010101 

1.1723521 

0.23000 

1.11217412 

115.117124 

715.17331 

111.3707 

722.110 

3.440 

70.44711 

1.1000100 

1.1723521 

0.21000 

0.27744431 

115.117124 

731.10517 

MO. 71 11 

147,311 

3.441 

70.44711 

0.1010100 

1.1723521 

0.21200 

1.32111103 

115.117124 

177.01404 

HO. 3012 

107.707 

3.440 

70,44711 

1. 00^0001 

1.1723521 

0.21100 

0.31143137 

115.117124 

117.14351 

137.1712 

711.141 

3.440 

70.44711 

1.0010100 

1.1723521 

0.21400 

1.30343741 

115.117124 

753.10333 

131.2407 

447.741 

3.440 

70.44711 

0.1100000 

1.1723521 

0.21000 

0.13773773 

113.117124 

477.74033 

134.3114 

511.H33 

3.440 

71.44711 

0.0000000 

1.1723521 

0.25000 

1.33271041 

113.117124 

373.31311 

132.3372 

173.110 

3.440 

70.44711 

0.0000000 

1.1723521 

0.25200 

0.74110031 

115.117124 

172.11277 

121.3112 

313.437 

3.440 

71.44711 

1.1000001 

1.1723521 

0.25300 

1.10000000 

113.117124 

147.25770 

I0I.724S 

221.371 

3.440 

70.44711 

1.1011000 

1.1723521 
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3 

DIAMETE1> 
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5 
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4 
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0.0017347173 
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7 

DIAMETER* 
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1 

DIAMETER* 
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IXTERATION> 

7 

DIAMETER* 
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IKTERATIOK* 

10 

DIAMETER* 

0.0011417731 

IKTERATION« 

11 
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12 
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13 
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0.0011241431 
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(• THIS PIOCRAH OETERHINES CERTIIN PIRAHETERS OP TriE EPFlUl 
(■ lESULTINC FROM THE PARASITIC HEAT LOAD GH A HELIUM DEVAR IN 
(• SPACE. THE SHIELD SYSTEM CONSISTS OF CONSTANT TEMPERATURE 
(• SHIELDS WITH A GIVEN DIAMETER AND PIPE LENGTHS. THE MASS FLOV 
(■ AND SHIELD HEAT EXCHANGE ARE THEN CALCULATED ALONG VITH THE FLOW 
(• VARIABLES. 


LAIEL II; 

CONST 

PUl.HlSmSS; 

l>XI77;(*J/KC-E*I 

I>l.i7; 

CP>S2IB; 

TYPE 

DIVISIOKS>ASRArtS..20Q01 OP REAL; 
'COVKr>ARRArCf..2G] OF REAL; 
ICOUKrsARRATC0..2l)l OF INTEGER; 
RCA10UT*FILE OF HEAL; 

IKIOUTsPILE or INTEGER; 


NUMVAlLS,NUNDIV,X,MrEST,PEEV: INTEGER; 

I . LOOPCOUNT. TAB , P IPECNT, VALLCNT, SPACECNT : INTEGER; 
EPPOUTM,POUT:TE2T; 

EPPINMi.REALOUT; 

EFPINH2: INTOOT; 

P.T.BHOrOIVISIONIi; 

I,Q2.H.PO,TO,Q,RE,F,NU,MDCT:OIV1SION5; 

erP,TV,lSPACC.LriFE,LVALL,aSHIELO:CDUXI; 

l.LTOTAL.OIA,LOVSOVNO,HIEOUND;REAL; 

TBSAVE .OSAVE .KDOTl , HDOTMAX .MAZEFF : REAL; 


Ri.TOf.Miji,Pi.iHoi.rt.ici,Qi,aii.Nvi:imi 


PlOCtOUlC IKDATA( IXDCZ : IKTCKl) ; 

ICCIN 

ir INDEX > I THEN 
ICGIN 

(• lEAD IN DATA *) 

VIITEIKC INPUT DATA AS PIORPTEO IT TEININALM; 

UIITEINC *); 

VIITELNI 'WEN ASEEO FOR A SERIES OF DATA SUCH AS THEM; 
niTElNCSHIElO TEKFERATUXES, INPUT THE SATA IN'); 
niTElNI 'ORDER FION TNE DEVAI TO THE OUTER lOUNDAIT*); 
VIITEINC '); 

VIITEINI 'ENTER THE NUNIEI OF UOIATION SHIELDS'); 

lEAOlN; 

lEAO(NVHVALlS); 

VEITEIEFFIloa.NUNVALES); 
niTELNCNURBER OF SHIELDS*' ,NUHVAlLS:S) ; 

VIITELNC ')' 

niTELNI 'ENTER NUNIEI OF SECTIONS FOI SOLUTION'); 

lEAOLN; 

lEAO(NUHDIV); 

VllTEIEFFINNZ.NUKDir); 

niTELNCKUKIEI OF SOLUTION SECTIONS*' ,NUKOIV:S); 
VIITELNC '); 

VXITELNI 'ENTER OIAMETEI (H)'); 

lEAOLN: 

lEAO(OIA); 

VIXTEIEFFINNI.OIA); 

VIITELNI 'DIAMETER (H)*',0IA:1S:F); 

VIITELNC '); 

VIITELNC INPUT DISTANCES lETVEEN SHIELDS (N)'); 

FOI I:*l TO NUKVALLSil DO 
ICCIN 
lEADLN; 

lEAO(LSPACEm); 

niTEICFFINNl.LSPACEin); 

END; {‘FORM 

VIITELNCCNTCI PIPE LENGTHS ALONG SHIELDS <M)'): 

FOI I: >1 TO KVXVALLS DO 
lECIN 
lEAOLN; 

lEAO(lVALLm); 

VIITEIEFFINMl.LVALLm); 

END;<*F01*) 

VIITEINCENTER TEMPERATURES OF RADIATION SHIELDS (X)*); 
FOI I:*l TO KUMVALLS DO 
lEGIN 
lEADLI; 
lEAO(TVm); 

VIITE(EFFINNl.TVm); 

END;(t FOI <> 

VIITELNCEN7ER MAXIMUM SHIELD EXCHANGER EFFECTIVENESS'); 
lEAOLN; 

ICAO(NAXEFF); 

VRITECEFFINMI.MAXEFF); 

VRITELNCMAIIMUN SHIELD EFFECTIVENESS*' ,MAXEFF;lS:f); 
VIITELNC '); 

VCITELNCENTER INITIAL STAGNATION T (X)'); 

READLN; 

REAOmm), 

VRITCIEFFINXl.TOni): 

VIITELNC INITIAL STAGNATION TEMPERATURE (X)*' ,Tim:I3;P); 
VHIIELNC '); 

7V(D]:«T0m; 

VIITELNI 'INTER INITIAL STAGNATION P (PA)'); 

RCAOLN; 

REAO(POIIl); 

VIITEIEFFIXNl.Fflm); 

VRITELNC INITIAL STAGNATION flESSUlE (PA)*',POin:lS:f); 
VRITEINC '): 

VRITELNC ENTER KISS FLOV APFROZIMATION (XC/S)'); 
lEADLN; 

RElOIMDOTm); 

VRITEIEFFINMl.MDOTCll); 

VRITEINCHASS FLOV (XC/S)*'.MDOTCn:lS:f); 

VRITELNC ')■ 

VRITELNC 'ENTER NUMIEI OF STEPS IN RUKCE-XUTTA METHOD'); 


ORIGINAL PAGE 
OF POOR QUALITY 


ICIDUI; 

ICIO(K): 

vimdrriKHi.M); 

niTtiKi'MUKiEi or snrs in iuNGMvnA<*.N:S): 
tNOIMF*) 

nsc 

(• Ku rill *) 

iceiN 

itsmirriNHi): 

Ksniimwa); 

ICiOdFriKNI.KimVUlS); 

lEiodrrnoii.KUKDiv): 

lUOdrriNMi.oii): 

roi I:>l TO KUHViLlS«l 00 REAOdmiOll.LSPACEIlDi 
roi I:>1 TO NUHVALLS OO BEAOdFFINKl.LVALim); 
rot I:>l TO KimVAlLS 00 lEAOdFFINni.TVm); 
lEAD(EFFINKt.HAZEFF); 
lEAOIEFFimi.TOm): 

TVm:>T0m: 

lEAOIEFFIKKI.POm); 

lEAOiEFFINMl.HOOTm); 

lEADIEFFINKZ.N): 

EN0(*CLSE<) 

SND;(« UOATA M 


PlOCEDtflE OVTDATA; 
lECIN 

VliTElKIEFFOimf.'LANFlOV 0VT?UT'.‘<l); 

WITEIHIEFFOVTH,' 

VIITElNIEFFOinN, ' '); 

WIITEIHIEFFOUTM,' *>: * 

VIlTElNdFFOUTH. 'VALUES OF IKFOTMd); 

WIITEINIEFFOUTH,' '); 

niTElNIEFFOOTK, 'NUNBER OF SECTIONS «' ;SO,NVNOtV;S); 

VIITELNIEFFOUTN.'DEVAR STACKATIOK TEKPEBATURE (E)>';S«,T0m:2l); 

VRITElNIEFFOVTH.'OEVAfl STAGNATION PRESSURE (PAU' ;SO,POC11:20); 

VIITEUKEFrOVTM.'NASS FIOV AFPIOZINATION dC/S>>*:SO,NOOTm:2f ); 

WiTElNIEFFOUTN, • ' 

VRITELNIEFFOUTH, ' 

VtlTELNCEFFOirnt, 'SHIELD nKFElATVRES (K) 

FOR I:sl TO KUHVALLS DO VRITELNIEFFOUTN. 'SHIELD TEMPERATURE' :4f, 1:2. :l,TV[n:iS;P)i 
VIITELNISFFOUTK, ' '}; 

VIITELNIEFFOOTM. ' 

VIITElKIEFFOirnt, ' '); 

VIITELN(EFFQUTN, 'VALUES OF OUTPUT* :I0); 

VRITELNIEFFOUTH, ' 

VIITELNIEFFOUTM. 'HASS FLOW (EC/S)-' : SB, NDOTIIOOFCOUNTI: 13:2); 

VIITELNICFFOUTH. 'TOTAL HEAT TRANSFER (V)-' :SO,aZ[NUKDIV];lS:F): 

VRITEINCEFFOUTH, ' '); 

VIITEINIEFFOUTX, 'SHIELD HEAT TRANSFER (V)';d)i 

FOR I;«l TO KUHVALLS DO VRITEINIEFFOUTN, 'SHIELD HEAT (V}':40,I:2,'-';l,QSRIElOni:lS:»): 
VIITELH(EFFOUTH,' '); 

VRITELNIErrOUTN, 'SHIELD HEAT EZCHANCEI EFFECTIVENESS' ;iO); 

FOR I:«i TO KUHVALLS DO VIITELN(EFFOUTH, 'SHIELD EFFECTIVENESS' :d, 1:2, '■' :l,EFFin:!S :t) 
VIITELN(EFFOUTH, ' '); 

VIITELHIEFFOUTH, ' 

VRITELNCFOR AN OUTPUT OF SECTION FLOW VARIABLES, INPUT 1. OTHERWISE INPUT I.'); 

READLN: 

REAOIPKET); 

IF PEET « 1 THEN 
BEGIN 

WRITEILFFOUTM, 'LENGTH':!#, 'MACH HUMBER': 13, 'STAC TEMP' : 13. 'STAC PIES':13); 
WRlTEdFFOUTH, 'TEMP' ;lfl, 'PRESSURE' ;10,'NUSSELT #' :10, 'RETHOLD f :12); 
VRITELNIEFFOUTM, 'SECTION HEAT' : IS, 'TOTAL HEAT' : IS): 

WRITElKIEFFOUm, ' '); 
rot I:>1 TO NUMSIY OO 
BEGIN 

VRITE(EFFOUm,' ' :?„ZtI] :i:3, ' ':5,Min:10:l,' ':5,T0m:ll:n; 

VRITEIEFFOUTK,' ' ;S,P#tIl: i#:S, ' ' :2,TII1:I:4,' ' :2,Pm;l:2); 

VRITEIEFFOUTM.' ' :2,XU[I1:I:3. ' ':3.BE(I1;10:S,' ':S, 0111:10:7); 
VRITEIMCEFFOUTM, ' ' :S,azm: 11:7) ; 

CND(*FOI<) 

CND:(<1F<) 

VRITELNC '); 

VRITELNCFOR A PLOT FILE OfTER 1. OTMSRVISE ENTER #'); 

READLN: READ(TAB); 

IF TAB > I THEN 
BEGIN 

VRITELNCFOR PLOT OF STAGNATION TLMFERATVRE ENTER 1.'): 
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MimMCFOl PLOT or STACN/.T10K PSCSSURC EKTCI 2.M; 
niniNC’roi plot or heat TiAitsrEi entei j.'): 

ItAOUilElOCTAI); 
roi 1:*1 TO KUKDIV DO 
lEClN 

CUE Tu or 

t: VlimH(POUT.TI[in; 

1: VllTCLM{POOT,PdEin; 

1: VIITCLN(POUT,Qim)i 

CMD(*CASE*) 

END(*rOIt) 

END;(<ir*) 

vumNcroi ixtciation intomation, ihpot i. otheivise ikput e/)i 

lEAOLN; 
lEU(PEn): 
ir PIET ■ 1 THEM 

roi I:«l TO loopcouht do 

VRITELKi EFFOinM. ' IHTERATION< ' ; 21 , I : ^ . 'HASS FL0V> ' ; 21 .RDOTE 11:23); 
VUTEINCOOTPUT STOIEO IM FILE EFFOUTN'),' 

DO:(* OUTDATA *) 


rWCTIOH PVI(N.P;IE1L):IEAL; 
(* SOLVES roi K*>P *) 

VU A.lEAl; 

IICIM 

l:>P*tMfH); 

m;.EXP(A): 

EMO:(* PWl i) 


FUNCTION f(H:IEAL];BEAL; ' 

(• SOLVES OFTEN OCCUMIKC FUNCTION *) 
IE62H 

l:>l.l«<I-l.l)*SQI{N)/2.l; 

INOi(> I *) 


rUNCTION TrNP(TO,H:IEAL):ICAL; 

(• CALCULATES LOCAL TEMPEUTUIE «) 
SECIN 

TDIP:>TI/l(Hi; 

END:(< TEMP >1 


FUNCTION riCS(P0,N:8EAL);IEAL: 

(« CALCULATES LOCAL PRESSUIE <) 

m A.C:1EAL; 

lEClN 

C:>I(H): 

1:>E/(1*E); 

PiES:>PI<FVl{C.A}; 

ENB;(< PIES t) 


FUNCTION OENSC{P.T:REAL) :BEAL; 

(• CALCULAnS LOCAL OENSITT M 
lEClN 

OENSE:>PMl*T); 

ENO;(i DENSE M 


FUNCTION MU(T:IEAL);REAL: 
m A,a,C.O:lEAL; 

IEC1H(<NUM 

A:il.l041t2A3S31E-0<i 

l:B|.l(Si77SM24IE-l7*T; 

C:>-3.232722F210i2E-ll*SQI(T)', 

0:b4.72S700273F22E-13<T<SQ1{T) ; 

HU.'iAtBtCiO; 

END;(>HU«) 


FUNCTION C0ND(T:1CAU:REAL; 

VAI A,B,C.O:lEALi 
lECINItCONSUCTIVm CALCULATION*) 
A:*I.I1S7213S71; 
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C:B-i.»H07Sm-l7*SQR(T)i 

0;>4.3ZII3I7(MI*T*SQ1(T): 

C 0 ND;a«i«C 40 : 

nOi<*coio)«) 


mCTIOM FIUIDTl(T:IEiL):ICll| 

111 i.l.C.O.E,r:ini: 

IEC1N(«PIAKDTI KUKBEI CilCUUTIOMM 

A;>l. 371237771; 

l:.MI33ini71«T: 

C:>>7.8I324332E-I3>SQ1(T); 

D;a4.3272a3BlE-l7iT<SQX(T); 

E:*-1.1I«1443E.17*SQ1(T)(SQ1(T); 

r:>l.l74l3l42E-t2*TtSQl(T;<SQS(T); 

FIUa>n;>i*l*Ca*EtF: 

END:(*FliKDTl*) 


rUNCTIOX HriIC(Hl.lU.Oir.l:IEU}:IEU; 

Vlt HFI: 1 EU; 

(« ISEO IN CASE 2 OF MCVTON <) 

IE61N 

RFI:>(Sai<H 2 )-SQI{Hn)/<K<SQI(Hl«H 2 )}«(E«l.l)/( 2 . 1 *C}<lN(l(R 2 )*SQI(Kl)/(l(Hl}*SQl(H 2 )))- 4 .l*F*L/D 

RF&1C:>NFI; 

EIOi(* KFllC M 


FVXniON DNF1IC(H1,H2:SEAI):IEAI; 

(t OSED IN CASE 2 OF MCVTON <) 

7AI &KF:ICAt; 
lEClN 

0HF:>2.l/(CtH2tSQl(X2))4(t«nyE*H2*Sai(l(Hl}}Mfn(Nl,4><FVI(l(H2).3)); 

0HFtIC:>0HF: 

IF ONF < 1.274 THEN (•FUNCTION UNSTAILE, lUCH*t IFFIOACHED*) ■- 
DHFIlC:*i.l: 

ED:(* DHFIIC •) 


FUNCTION IfITSST (H . 0 , ?aOT , F I , TO : READ : HEAL : 

(• USED IN CASE 3 OF HEVTON •) 

VU A,C:1CAI; 

ICCIN 

A:bI(R): 

C:><I*l.l)f(2.|t(t.l-l}): 

RFRST:*{R«FO*Ft<SaR(O)M.0)tSaiT(K/(E*TintFyR(A,C}-RDOT; 
END; (I RFIST •) 


FUNCTION 0HFKST(N. 0 . FS , T8 : HEAL ) :HEAt ; 

(• USED IN CASE 2 OF XEVTOK *) 

m A,C,C:IEAL; 

lEClN 

A:>I(R); 

C:«(K«2.l)/(2.|i(l.l-C)); 

E:«(3.|tE-t. 3)1(2. 0*(1.6-E)); 

DHF15T:«(FI>SQR(0)/4.3)«F0tSaRT(E/(R*TS))<(FVK(A.C)-SQI(R)t(K«l.l)y2.g<FVR(A.E)); 
CXD(* OKFRST •); 


FUNCTION KEVTONCEETrlNTECEX; Ml,O,F,l,lTOT,F0.TI:lCAL):HEAl; 

7U KKCV,KOLO,OIFF,Ol.02,A.MDOT.OKF:HCAL; 

COUNT: INTECEI; 

IECIK(* KEVTOH *) 

DIFF:>I.I: 

CASE CET OF 
1: lECIN 

(• CALCULATES EXIT RACE FOR SECTION VI7R FRICTION ONLT •) 
X0LD:<I.S1*H1; 

COUKTrit: 

VHILE (OIFF >« 4.30003) AND (COUNT ( IS) 00 
lECIN 

DKF;*OKFRIC(R1,NOIO); 

KNEV:aOLO-NriIC(Hl.HOLO,0,F.L)/ORF; 

DIFT.-aBSlKOLO-HNCV); 

IF ONF « 1.0 THEN ('FUNCTION UNSTARLC. EXIT LOOF*) 
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IICIN 

Dtrr:«MOI01; 

IDKV:il.l; 

E0i(MF«) 

comr.-aCOUKr*!; 

MtO:*KKEV: 

C0(« VKIII M 
DID:(*CiSE 1*1 
1: lEClK 

(• CUCDUTES INLET RICE K/KEEI M 
ND0T:>L: 

WILE (Oiirr )■ i.Mflit) 00 
lECIN 

HXEV: >NOlO-MriST(KOLO , 0 . KDOT . r 9 , Tl ) /DMriSTimLO . 0 , ro , 70 ) : 
DirF;>ilS(HQlO-HNEV); 

K0L0:>HNEV; 

END(« VRllE M 
EXDI« 2 •) 

EKD;<* CISE •) 

NEVrON:>NOlD; 

END;(* NEVrOX t> 


nOCESUIE INlETKlCR<O.RDirr,FO,T 9 :lCAL: YU N.T.F.IXOilEU): 
(a CUCVIATES INLET CONDITIONS <) 

TAI AJEALi 
lECINI* IKLETMCX «) 

H; -NEVrONd ,N, 0. 9 .ROOT. 9 J . Ft .Tl i : 

T;»TENF(TO,N); 

P:iFIES(PI.N>; 

UO:>OENSE(PJ); 

IXD;<* INLETXACH *) 


PIOCEOUIE FHICTIOH(TO,OIA,F,1:IEAL;VU NJ,P.lHO,PO,a:lCAL); 
(* FOE FLOV VI TH FRICTION ONLT «) 

YAI A.C.REVH:SEAL; 
lEClH 

NEVH:-NEVrON( 1 .N.OIA.F ,L .LTOTAL.I , I) ; 

A:«I(NEVM)/B(N); 

C:>(X«1.I)/(2.I>(E>1.I)); 

PI:sPO*NmCVNtpVR(A,C): 

R:«XEVH; 

T:»nNP<TO,NJ; 

P:«P1C5(P9,HI; 

IRO:>DENSE(P.T)i 

Q:.l; 

EXS;(<FIICTION«) 


FUNCTION STACTEKF(TO,TV,NU,IF.XDOT,L:REAL);lCil; 

(t CALCULATES NEV STACHATIOH TEKPERATVXE FOR FLOV VITH HEAT TRINSIEI «) 
lEClN 

STACTENP : xTVt ( (TB-IV) / EZP < ( NU*P I a EF* L ) / (NDOT*CP ) ) ) ; 

END:(<STACTENP<) 

FUNCTION IFUNC(X,T9,1V,FR,NU:lEiL; KE7; IKTECCI) :REAL; 

(• SOLVES AN OFTEN USED FUNCTION IN FEOCEDURE HEATTIANS 

TAI A.C,D,E:1EAL: 

lECIN 

AixRtKR) /(2.0<( 1 . S-SQR(N) ) ) ; 

C:«(l.l«E*Sai(NmT9; 

D:sE*SQI(R); 

IF EET • 1 THEN 
lECIN 

E:«2.i/nV-T0); 

ENDdlFt) 

ELSE 

lECIN 

E:*UJ*FII(NU«(7V-Tf))i 

END:(tELSE*> 

IFUNC:iA<(OtEfC); 

EXD:(< IFUNC •) 


riocnuic licmiiMscTv.i.o.norjcu; nciktcch: m H.Ti,ri.T«MHO,».r,o.ax,xu;im>: 
i* aumm riov vAiuaiis foi riov vitx fiictiok 110 meat hansfei <) 

VAl 

Cl ,C2 ,C1 ,C4 .snp , A .C ,NEVl.FAC70l ;ICAl i 
Tl I ,TI2 .Ml , HI . PMF , OCIR: IE Al ; 

[ir,I,nSTt!KTECUi 

iieiM 

nx:>i: 

nST:*l; 

<t FIAC lEI IF FLOW IS rJIBULENT t) 

It m ) t.H THEN IET:il; 

(• eUCVlATE NEW STACNATION TEHFEIATUIC M 

Pl:*PlAKDTl(T)i 

ir:>CONO(T); 

TI1:-TI; 

TIE: >STAGtEKP(TO I ,TV, XU , IF ,;tDOT. 1) ; 

!P (1V-T82> < 0.0001 THEN 

(* VALL TEKPEIATURE APPIOACIEO, EXIT PIOCIAN *) 
lIClM 

A:i-(n»T*CPmPI>NU'IF); 

NEVl:«A<lN(l.02); 

rA7rOI:il/NEVl; 

niTElNCTHE FLUID TEKPEIATURE AFPIOACHES THE VALL TEKPEIATUREM; 

VIITEINCAT AN IKOETERKINATE LENGTH. INCEEASE THE NVNBCI OF'}; 
niTElNl 'SOLUTION SECTIONS. SUGGESTED INCREASE IT FACTOR OF'. FACTOIDS, 
VIITELNCOI MORE.'); 

G01D IS: 

END;<<IF*I 

’ (« STEPWISE SOLirriOM USING lUNGE-EUTTE METHOD M 
I* eUCULATE STS:P SIZE <1 
STZP;>AIS((TI2-TIU/N>; 

I:>1: 

VIILE (I <> N) AND (TEST • 0) DO 
BEGIN 

Cl : iRFUNC (H.TO , TV. FI ,HU, RET ) ; 

IP Cl < 0 THEN Cl:«0;(tFUKCTION UNSTABLE*) 
C2:.|FUN(i((H*Cl»STE?/2.0>,(TB4STEP/2.0),TV,FH,MU,RCT): ' 

IP C2 < 0 THEN C2:il;(*FUNCT10N UNSTABLE*) - 
CS:ilFUNC((K4CZ*STEPI2.0),(T0«STEr(2.0).TV.?R.NU,XET): 

IP C3 ( 0 THEN C3;iO;(*FUNCT10N UNSTABLE*) 
C0:iirUHC(H*C3*ST£P,TB4STE?,TV,PI,XU,XEI): 

IP CO < 0 THEN Cl ;>0;( (FUNCTION UNSTABLE*) 

OElH;iSTEPU.0*(Cl4l.0*C24l.C:C34CI); 

R:«K40ELH; 

TI:>Ti4STEP: 

I:>UI: 

IP K ) 1 THEN ((FUNCTION UNSTABLE, EXIT LOOP*) 

BEGIN 

TEST;.I; 

l;>l.l; 

END(*IF*) 

ED:(* WHILE *) 

IP (Cl > I) OR (C2 > 0) OR (C3 * 0) OR (CO * I) THEN 
((FUNCTION UNSTABLE. EXIT LOOF*) 

BEGIN 

TEST:«1: 

H:ii; 

XNO;(*IF*) 

RS:iK; 

A:>l(N2)i'B(Hl); 

C:>(ltl.l)/(2.9*(X-l.f))-, 

PI;«PB*(Kl/K2}(SQlT(T02/T01)iFWRa.C); 

P:>FRCS(F0,H2}; 

T:.TEHP(T02.H2>; 

IRO:*OENSE(P.T); 

IE:sO.I*NDQT/(FI*D*HU(T)); 

IP IE (« 2000 THEN 

BEGIN (( LAHIKAR FLOW *) 

F:itl.0/RE: 

NU:s3.1<: 

E)a(*IF() 

ELSE 

BEGIN (* TURBULENT FLOW *) 

F;>I.III/FV1(RC,0.2): 

NV:sO.013*FVR(RC.O.I)>PVI(PR.I.I); 

END:(*ELSE*) 

Q:iCP*KDOT»(T01.Tfll); 

QX:«QZ4a: 


O0i(« ukTmta •> 
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nocnuic SECTIOMnWIlLiIKTECEl; VEI HSITilKnCER; LS,lV;COUKT; VIR X:DITISIOMS); 
(« THIS PROCESDRE DIVIDES THE PIPE LENGTHS INTO MOHOVERUPPINC SECTIONS *) 

m 

IESl,LTDT,L:REAl; 

cooKri .comn , i , j.ietktecei; 

IRC» (* SECTION i) 

(• CRICVUTE KOHIKiL SECTISN SIEE •> 
lTOT:>li 

roi l:>l TO NVIll DO lTOT;>LTQT«LV(II: 
roi t:>l TO KVULtI DO LTOT:>LTOTUSm; 
l;-lTtrr/KUM01V; 

(• CAECUIATE PIPE SECTIONS M 
LTOT:>l; 

Zm:>l; 

J:-l; 

EtT:>l; 

COUXTt:«l; 

COUNT!; at; 

roi i:>i TO 00 

lECIN 

III;«IEf«l; 

It EZr ■ 1 THEN 
lECIN - 

lTOT;>lTOT{lStCOUNTi]; 

COaNTl:«CDUKTUl: 

END(«ir*) 

ELSE 

lECIN 

LTOT;<LTOTflVCCOUNT!l: 

eOinnn:iCOUNT 2 *i; 

IET:>0: 

EN0;(<EISE*) 

VRllE XIJ1 ( LTOT OO 

lECIN *• 

XIJl:aX[J-ll*L; 

END;(<VHIlEt) 
ir XXJ] > ITOT THEN 
lECIN 

IESL:-ITOT-X[J-1I: 

IU]:aX(j.ii*SESL: 

ENO;(iirt) 

CND;(<fORt) 

NDir;>J; 

IINDmi];ilTOTfl; 

END; < (SECTION*) 


PROCEDVRE SETBOVKDdETrlHTCCEl; ?U HII,COVI.HDQT:ICIl): 
lECINItSCTIOUNOi) 
ir EET ■ I THEN 
lECIN 

(OVI:>M!IOT; 

END(((F() 

ELSE 

lECIH 

III:>KDOT; 

EHD;<(EISE*) 

NDOT:><2.t*lOVIiRII}/3.l: 

END:((SSTBOUND*> 


PROCECUIE INHEASSICH; 

((DEFINES VALUES FOR USE IN PAIANETER LIST *) 
lEClNI* INREASS ICK() 

L:>im<XIMJ; 

Nl:>Hm; 

Tll:sTICIl; 

PllisPOCIIi 

Tli.TlII; 

Pi:>P[Il; 

RH01:*IHOm: 


n:irni: 

ici:>ir;n; 

ni:>iW(ni 

Qitiom; 

OlluQim; 

ED;(*inElSSICM» 
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riOCCOVIC OtmCASSICK; 

(• DCriMSS NIV VALUES OITAINCO FlOH rAJtANCni LIST •) 
ICClK(*OUnEASSICK«) 

l:«Uli * 

Tltn:cTll: 

Tm:*Tl; 

Uom;>IH01: 

rm:>rii 

am:iQii 

QHn:«QXl: 

ao;(*ounEAssicM«) 


Piocnuit INITIAUXt; 

(• PKtrUES VAIIAIIES FOI USE IN MAIN PlOdJUl <) 

' TU 

jMncn; 

A,C:IEAt; 

IECIK(«INITIALIEEM 

KDOrm-iO; 

(fOEtEIKlKE NAXIKUH KISS riQV«t 
A:>LltMISMI-La); 

C:><I»M)I<2.0M1.M)); 

HD(mUX;>(PI*SQX(OIA)/1.B)<POtl]tSQBTtI/<l*TO(nn<PVI(A.C): . 
{•CXECX GIVEN NASS FIOV AND BEOUCE IF NECESSARY*) 

IF KOOTII] ) ItDOTKAX THEN KOOTIl];>«.7*KOOTXAX; 

(« CAtCUUTE TOTAL FIFE LENGTH ») 
lTOTAL:>l; 

FOI I;il TO KURVALLS DO LTQTU:>LTOTAL«LVALLm : 

FOI I: tl TO KUKVALLStI 00 LTOTAL:>LTOTALUSFACEin; 

(«CAICULATE SECTION LENGTHS*) 

(■INITIALIZE SECTION ARIAT*) 

SECT! ONINVKVAILS ,NUKD I V , LSPACE , LVAll , X ) i 
(*STA1T LOOP COVNTEl*) 

IOOFCOUKT;>B: 

(* INITIAL NACH NUKBEI GUESS*) 

NdlriO.ll; 

(* INITIALIZE END KACH NUKBEI SO THAT IT KAY BE USED FOI TEST*) 
NINUKDIVI;>O.S; 

(*INIT1AFJZE BOUNDS *) 

HIIOOND:*HDOTKAX; 

IOVIOVNS;«M; 

END:(*(NITIAIIZE*) 


FIOCXOUIE lElNITIAL; 

<* PIEPAIES VAIIAILES FOI USE IN KAIN PIOGIAI *) 
TAR JzIKTECEl: 

IEC1H(*REINITIAI*) 

(* INITIALIZE COUNTEIS*) 

I:sl; 

FIPECirr:>(; 

VALic)rr:«i: 

SFACECNTzati 

{*INCIEMEKT LOOP COUNTER*) 
lOOFCOUNTziLOOFCOUHT^li 
(*CALCULATE ENTRANCE CONDITIONS*) 
LFIFEmzsLSFAGCni; 

QXm:>l;Qm:«0; 

lEtn ;>4. 0<KDOT(lOOFCOUNnMPI*Kn(TS(ll)*OIA) ; 

IF icm <■ loog THEN 

BEGIN (t lAKINAE FLOV *) 

Fm:<UJ/RE(n; 

NVm:>J.di 

END(*IF*) 

USE 
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ICGIH (* TUHUIEKT FIOV *) 

r(ii:>M«/rvi(icnM.2): 

IUtU;«M2]tpwR(icnM.I)«PVI(PIIUaTl(T0ni).M); 

CN0;(ICISE*> 

(•CllCVUTK INICT COKDII7IONS*> 

iximucK(ou,Hi)OT(ioopcotmTi,pon]ji(n.RCi] jni.P(ii.iHoni)! 
(tlMITIUIZE comoi CAM*) 

HnSTl.1} 

aiP:(*iEiNiTiu*) 


lteiX(« MIN PIOCIIN M 

UiTElHCTO INTEI OATi, EKTEM. OTHElVISE.EKTEt l.‘); 

lEAOLNi 

lEAOiTAII: 

IIOATA(TAI); 

INITIALIZE: 

lEPEAT (* UNTIL THE FLOV IS CHOIED *) 

ICINITIAL: 

VIILI (I < KVKOIV) AND (HCII <« 1 . 17 ) 00 
lECIN 

(• lECIH FRICTION FLOV CALCULATIONS •) 

VIILt (ZLII <> LFIPEIPIPECNT]) AND (NIII <• 4 .F 7 ) DO 

mix 

INIEASSICN: 

FIICTION(T 01 . 0 IA,F 1 .L,H 1 .T 1 .P 1 ,IH 01 .POI.Q 1 ): 

OUnEASSlCN; 

CND:(*VHILE*) 

IT VALICNT (■ HUKVALLS THEN 
lECIX 

PIPECXT:«PIFECNT 4 i: 

SPACECNT:>SPACECKT«t; 

lPIPECPIPECNT]:>lPIPE(PIPECNT-n»LVALLIVALLCXn: 

END:(>ir*) 

(• lEGIN HEAT TRANSFER CALCULATIONS •) 

TOSAVEtiTItl]; 

EmVALLCKT]:«l; 

QSAVC:iQim; '• 

WILE (EFFLVALLCXn ( NAZCFF) AND (ZIII < LPIPEIPIPECKT]) AND (NCI] (■ 1 . 77 ) DO 
lEGIX 

INIEAS 5 ICK: 

HCATTRAHS(TVCVALLCNTl,L,OIA,)(DQnLOOPCOtfNTI.X,Hl.r 01 .P 01 .TI,Pl.lROI,REI.ri,ai,QZl.NUI) 

0 U 71 EASSICN; 

EFFIVALICNTI : > (TO ( 1 1 -TDSAVE) / (TVIVALLCNTI-TSSAVE) ; 

END;(tVHILEt) 

OSIIELDCVALLCNT) : < 0 Z C 1 1 -OSATE ; 

IF VALICNT (« NUNVALLS THEN 
lEGIN 

PIPECXT:>PIFECNTtl: 

VALLCNT:>VALLCNTtI: 

LPIPE(PIPECNT]:>lPIPELPIPECXr-U«LSPACELSPACECNT); 

FNOCMF*) 

END:(*VHILE«) 

IF NUXDIV ) I THEN NTEST:*!; 

HDOTI;iKDOnLOOPCOUKTl; 

SETIOUND(KTEST,HIBOUNO,LOVSOtJNO,n)OTl) ; 

HDOTL LOOPCOUKTf n : «ND0T1 ; 

UNTIL (KINUIDIV] ) 1 . 77 ) AND (NLNUKDm ( 1 . 03 ) 01 (LOOPCOUNT ■ ID; 

OUTDATl: 

If: EKD.(* KAIN FROCBAH >) 
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nw. 

TO ENTER DATA, ENTER 1.' OTHERWISE* ENTER 0. 

» I 

INPUT DATA AS PROMPTED BY TERMINAL 

WHEN ASKED FOR A SERIES OF DATA SUCH AS THE 
SHIELD TEMPERATURES, INPUT THE DATA IN 
ORDER FROM THE DEWAR TO THE OUTER BOUNDARY 

ENTER THE NUMBER OF RADIATION SHIELDS 
? 9 

NUMBER OF SHIELDS- 9 

ENTER NUMBER OF SECTIONS FOR SOLUTION 
f 129 

NUMBER OF SOLUTION SECTIONS- 129 

ENTER DIAMETIH (H) 

? O.OOU 

DIAMETER (H)« O.OOIAOOO 

INPUT DISTANCES BETWEEN SHIELDS (M) 

* O.OIY 
? 0.021 
? 0.024 

* 0.024 
? 0.034 
f 0.039 

ENTER PIPE LENGTHS ALONG SHIELDS (M) 

? 0.032 
f. 0.01 
t 0.008 
? 0.004 
? 0.038 

ENTER TEMPERATURES OF RADIATION SHIELDS (X) 

? 21 .7 
I 37.7 
» 49.4 
? 101.1 
7 144.4 

ENTER MAXIMUM SHIELD EXCHANGER EFFECTIVENESS 
7 0.98 

MAXIMUM SHIELD EFFECTIVENESS- 0.780000000 

ENTER INITIAL STAGNATION T (K) 

7 4.23 

INITIAL STAGNATION TEMPERATURE <K)- 4.230000000 

ENTER INITIAL STAGNATION P (PA) 

7 2100.0 

INITIAL STAGNATION PRESSURE (PA)- 2100.000000000 

ENTER MASS FLOW APPROXIMATION (XG/S) 

7 0.00000143 

HASS FLOW (KG/S)- 0.000001430 

ENTER NUMBER OF STEFS IN RUNGE-KUTTA METHOD 

7 4 

NUMBER OF STEPS IN RUNGE-KUTTA- 4 

FOR AN OUTPUT OF SECTION FLOW VARIABLES, INPUT 1. OTHERWISE INPUT 0. 
» 1 

FOR A PLOT FILE ENTER 1, OTHERWISE ENTER 0 
7 0 

FOR INTERATION INFORMATION. INPUT 1. OTHERWISE INPUT 0. 

7 1 

OUTPUT STORED IN FILE EFFOUTM 

1.824 CP SECS, 101444B CM USED. 

/ 
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nK.crrouTH 


umov ovnvT 


38 


fiivcs or ixpoT 


KVMin or SCCTIOMS « 
OEVll STICKATION TCHPERATURC (E)> 
OEVAI STACKiTION PRESSURE (?i)> 
USS riOV iPPlOXIKiTION (1C/S)> 


121 

4.23fllOOO«SegE*OI9 

2.1000SttOOO0E«0O3 

i.i3i(ioooieooE-ai4 


t 


SHIELD 

SHIELD TEHPERATURE t> 
SHIELD TEHPEIATUIE U 
SHIELD TEMPERATURE 3> 
SKPJID TEMPERATURE 4> 
SHISLD TEMPEUTURE S> 


TEHFEUTURES (1) 
21.790000001 
37.700000000 
43.400000001 
101.100000000 
144.400000001 


UALDES or OUTPUT 


I 

I 


I 


HASS riOV (KC/S)- O.OOOOOOISI 

TOTAL HEAT HANSFEI <V)« 0.4301S111S 

SHIELD HEAT TRANSFER (V) 
SHIELD HEAT (V) la 0.04971S14S 

SHIELD HEAT <V) 2a 0.042077032 

SHIELD HEAT (VI 3a 0.071471217 

SHIELD HEAT (VI 4a 0.073742177 

SHIELD HEAT (VI 3a 0.33 41 14403 

SHIELD HEAT EXCXAKCEI EFFECTIVENESS 
SHIELD EFFECTIVENESS la 0.043247702 

SHIELD EFFECTIVENESS 2a 0.443717702 

SHIELD EFFECTIVENESS 3a 0.434301374 

SHIELD EFFECTIVENESS 4a 0.300103174 

SHIELD EFFECTIVENESS 3a 0.704241774 


r 

LENCTH 

RACE HUMBER 

STAC TENT 

STAC PRES 

TEMP 

PRESSURE NUSSELT 1 

RETNOID 0 

SECTIOM HEAT 

TOTAL HEAT 

A. 


1.01473147 

4.230000 

2200.19010 

4.2277 

2077,417 

3.440 

427.31017 

0.0019000 

1. 1000000 


0.00403 

0.01470744 

4.230000 

2073.04477 

4.2277 

2072.443 

3.440 

427.30017 

1.0000000 

1.9000000 


0.00407 

1.11413401 

4.230000 

2014.11224 

4.2277 

2013.727 

3.440 

427.30017 

1.9000001 

l.iOIOOlO 

} 

0.00110 

1.01411414 

4.230090 

2077.20237 

4.2277 

2071.111 

3.441 

427.30117 

1.0000000 

0.0000000 

i 

0.01012 

0.01473342 

4.230000 

2072.31304 

4.2277 

2071.727 

3.440 

427.30037 

1.0090000 

I.IOOOOOO 


0.01214 

0.01471321 

4.230000 

2043.43022 

4.2277 

2043.043 

3.440 

427.30017 

0.0090900 

0.9000000 


0.01417 

0.01303311 

4.230000 

2031 .40711 

4.2277 

2031.217 


427.30017 

0.9000900 

1.1000090 

1 

0.11417 

0.01301307 

4.230000 

2031.70772 

4.2277 

2031.371 

3.’440 

427.30017 

1.0000000 

0.9000000 


0.91022 
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APPENDIX B: CALCULATION OF INLET MACH NUMBER 

The stagnation conditions of the dewar are equivalent to the entrance 
stagnation conditions. The continuity equation for one-dimensional flow 


states that 

ifj - pVA. (B.l) 

The definition of Mach number is 

M » V/C(yRT)^''^] (B.2) 

Combined with (B.l), this yields 

m « (yRT)^/^ pMA (B.3) 

The density may be expressed in terms of the stagnation density and the Mach 
number. 

P = P„n + C(y - l)/2] (B.4) 

The ideal gas equation states 

Combining Eqs. (B.4) and (B.5) yields 

P = [Po/(RTg)]D + (y - l)/2 (B.6) 

Similarly, T may be expressed in terms of M and Tq 

T = yi + (y - I)/2] (B.7) 

Combining Eqs. (B.3), (B.6), and (B.7) yields 

m = (t/RT^j) AMP^n + [(y - l)/2] (B.8) 


We wish to solve for the Mach number. Since the equation is nonlinear in M, 
the Newton-Raphson method of solution is invoked. Starting from some initial 
condition, new values of M are calculated iteratively until convergence is 
obtained. Using the expression 

«new = «old- (B.9) 

where 

f(M) = (y/RTq)^/^ AMP^n + C(y - l)/2] - ih (B.IO) 


- ORIGINAL PAGE IS 
OF POOR QUALITY 


42 


and 


f(M) . AP^(r/RT^)’^^ B<C-r+l]/2[(l-Yj)p _ m^C(t+1)/2] B’’} 
B is- an often used function of M and y, 

B - 1 + C(y - l)/2] 


i 

r 


\ 

i 


I 


I 


1 


(B.n) 

(B.12) 
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APPENDIX C: FRICTION FLOW MACH NUMBER CALCULATION 

The differential equation relating the change in Mach number to the 
change in distance is 

dM/M - (CyM^ B]/[2(1 - M^)]}[(4f dx)/D] (C.l) 

Solving with the imposed boundary conditions of M » M^ at x * and M » M 2 at 
X » X 2 yields 

[4f(x2 - x,)]/D » {(M^ - nf m|) 

+ C(t + r)/2r]*lnC(M^ B,)]} (C.2) 

Solve for M 2 with the Newton -Raphson method using 
f(Hj) = [(H^ - Mf)/(r nf M^)] 


+ C(t+1)/2yD ln[(H^ B2 )/(h|b,)] - 4fL/D 

(C.3) 

= [ 2 /(r» 2 ) + {[(-r+U B^)J 

(C.4) 


where L = X 2 - 
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APPENDIX D: STAGNATION TEMPERATURE CALCULATION FROM AN ENERGY BALANCE 

The heat transferred to a fluid as It travels through a pipe within 
isothermal walls is 

m dq * h(T^ - T^^) irD dx (0.1) 

The differential change in heat of an ideal compressible gas is a function of 
the stagnation temperature. 

(D.2) 

Combining Eqs. (D.1) and (0.2) yields 

^ Cp - 1^0 dx (D.3) 

The adiabatic wall temperature of a laminar flow is equal to the stagnation 
temperature. For turbulent flow, the slope of the temperature profile of a 
fully developed fluid is too steep to apply a recovery factor; therefore, the 
adiabatic wall temperature will be assumed to be also equal to the stagnation 
temperature for turbulent fluid. Equation (D.3) may then be expressed as 


^ dT„ = h(T^ - T^) If D dx (D.4) 

p 0 W O' 

One may recall that the Nusselt number is defined as 

Nu = hD/k^ (D.5) 

Equation (D.4) may be then expressed as m Cp dT^ = ir Nu k^(T^ - T^) dx. 

Integrating Eq. (D.6) with the boundary conditions 
Tq = Tq^ at X = x^ and T^ = T^g at x = Xg yields 

^02 “ “^w “ ^*^^w - ToiMexp(ir Nu L/ifi Cp)]} (D.7) 


The exit stagnation temperature of a length L of pipe is therefore a function 
of known quantitites. 
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APPENDIX E: HEAT TRANSFER WITH FRICTION 
The change in Mach number in a pipe is a function of the change in length 
and the change in stagnation temperature. 


a 

dM y M^B 4f dx . (1 + yM^)B 

(E.l) 

t 

From the energy balance in APPENDIX D, we have 


\ 

*> Cp dT„ . p V n/4 Cp dTp . n(T^^ - T^) » 0 dx 

(E.2) 

i 

i. 

One may express these relations in terms of dimensionless 

parameters. 

^ ' 

■ "^ 0 ^ = C4NUp/Pr(Rej3)] dx/D 

(E.3) 

i 

For laminar flow. 


■ 1 

Rejj - 16/f 

(E.4) 

i 

i 

[dTp/fTw - Tp)] - [4 %]/CPr(16/f)] dx/D 

(E.5) 

i 

1 

(4f dx/D) - (16 Pr/NUg)CdTp/(T^ - Tp)] 

(E.6) 

1 

i 

For laminar flow the differential equation is then. 


i . 

1 

dM .Y M^ B 16 Pr . . (1 + y M^) B 

” 2(1 - M^) NUjj ■ ^0 2(1 - M^) 

(E.7) 

1 

For turbulent flow, one may invoke the Reynolds analogy 


A 

St = f/2 

(E.3) 


and 

• 

1 


(E.9) 

i 

• 

^V^^w “ ”^ 0 ^ = dx/D 

(E.IO) 

• 

So that for turbulent flow. 

• 


dM Y B ^ ‘^’’"o (1 + Y M^) B ^'''o 

VT7TV 

(E.ll) 


iL 
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Equations (E. 7) and (E.ll) are solved numerically with the Runge-Kutta method, 
using about four Iterations, as shown in APPENDIX F. 
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APPENDIX F: RUNGE-KUTTA SOLUTION 
The basic Runge-Kutta equation is rather simple. 

Vl * "n * + 2C2 +2 Cj + c^) (F.I ) 

where 

'1 • 

<=2 ■ f(^on * "n ^ 

'3 ■ f(^on ^ \ * '2 

<=4 ■ ''(^on + ^^0* "n + =3 (''•2) 

The function F is simply a form of the differential equation we wish to solve 

for 

F(To„. «„) = dM/dTj„ (F.3) 

For laminar flow, 


c dM. 


MB . r JZ t. 16 Pr . N , 1 + V M.^t /,- 

- 77 ^ - Tj) * (f.4) 


For turbulent flow, 


dM. 


BM 


dT 


0 2(1 - M^) 


X <t-^) * 


(F.5) 
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APPENDIX G; FLUID PROPERTIES 

All curve fits shown below were performed by a standard IMSL Math Library 
routine on the University of Illinois Cyber 175, They are valid only In a 
temperature range of 4.2 K to 350 K. 

Thermal conductivity; 

» 0.0 157218071 + 6.3416525*10"V - 6.99407521*10“^^ 

-10 3 

+ 4.52003096*10 *T (G.l) 

Vi scosity ; 

U » 1.10419265531*10“® + 1 .160677594648*10“^*T 

- 3. 222922924062*1 0“^°*T^ + 4. 725700293922*1 0“^^*T^ (G.2) 

Prandtl number: 

Pr » 0.570232994 + 0. 0055872291 *T - 7. 58324532*1 0‘®*T^+ 4. 32720501 *10‘^*T^ 

- 1.1 1661445*1 0“®*T^ + 1 . 07403542 * 1 0 "^ 2 *j 5 


(G.3) 
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APPENDIX H: MACH NUMBER CALCULATION WITH. AREA CHANGE 
In most texts, the area ratio is given in terms of the area at the 
throat. A*. The following merely algebraically solves for some arbitrary 


A, /A* • n/M,) {[2/(t + 1)] (H.l) 

Aj/A* - n/Mj) (C2 /(y + 1)] (H.2) 

Aj/A, • (H.3) 

Solve for ijg with a Newton-Raphson method. 

f(M2> - CMj/Hj) - Aj/A, (H.4) 


r (M,) 


{Cy+1]/C2(y-1)]) ^ , B 


(H.6) 


